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ABSTRACT 

We present rotation curves for 19 early- type disk galaxies (SO - Sab). The galaxies 
span a B-band absolute magnitude range from —17.5 to —22, but the majority have a 
high luminosity with Mb < —20. Rotation velocities are measured from a combination 
of Hi velocity fields and long-slit optical emission line spectra along the major axis; 
the resulting rotation curves probe the gravitational potential on scales ranging from 
100 pc to 100 kpc. 

We find that the rotation curves generally rise rapidly in the central regions and often 
reach rotation velocities of 200 - 300 km/s within a few hundred parsecs of the centre. 
The detailed shape of the central rotation curves shows a clear dependence on the 
concentration of the stellar light distribution and the bulge-to-disk luminosity ratio: 
galaxies with highly concentrated stellar light distributions reach the maximum in 
their rotation curves at relatively smaller radii than galaxies with small bulges and a 
relatively diffuse light distribution. We interpret this as a strong indication that the 
dynamics in the central regions are dominated by the stellar mass. 
At intermediate radii, many rotation curves decline, with the asymptotic rotation ve- 
locity typically 10 - 20% lower than the maximum. The strength of the decline is 
correlated with the total luminosity of the galaxies, more luminous galaxies having on 
average more strongly declining rotation curves. At large radii, however, all declining 
rotation curves flatten out, indicating that substantial amounts of dark matter must 
be present in these galaxies too. 

A comparis o n of our rotation curves with the Universal Rotation Curve from 
I Persic et al.l (|l996l ) reveals large discrepancies between the observed and predicted 
rotation curves; we argue that rotation curves form a multi-parameter family which 
is too complex to describe with a simple formula depending on total luminosity only. 
In a number of galaxies from our sample, there is evidence for the presence of rapidly 
rotating gas in the inner few hundred parsecs from the centers. The inferred cen- 
tral masses and mass densities are too high to be explained by the observed stellar 
components and suggest the presence of supermassive black holes in these galaxies. 

Key words: galaxies: spiral - galaxies: lenticular - galaxies: structure - galaxies: 
fundamental parameters - galaxies: kinematics and dynamics - galaxies: haloes 



1 INTRODUCTION 

Rotation curves are the prime tool for studying the mass dis- 
tribution in disk galaxies. In normal, unperturbed galaxies, 
gas moves on circular orbits around the centre, so measure- 
ments of the circular velocity can be used to yield the en- 
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closed mass at different radii. The study of the shapes of ro- 
tation curves therefore gives important insight into the over- 
all distribution of mass in disk galaxies. Hi rotation curves 
in particular are useful, because they probe the mass dis- 
tribution to much larger radii than can be achieved with 
optical data and reach to the regions where dark matter 
dominates the gravitational poten tial. In fact, it was the dis- 
covery, first made in the 1970 's (|Rogstad fc Shostaklll972l : 
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iRoberts fc Whitehurstl 1 19751 ; iBosmal Il978l ; iBosmal Il98lh . 

that Hi rotation curves stay flat till the last measured 
points, well outside the optical disk, which gave the final, 
irrefutable evidence of the presence of large amounts of 
unseen matter in galaxies (Bosma 198ll; v an Albada et al.l 



Il985l ; lvan Albada fc Sancisilll986l ; lBegemanlll987i r 

A long standing question concerns the relation between 
the shape of rotation curves and other properties of indi- 
vidual galaxies. It has been known for a long time that 
the shape of a rotation curve is strongly coupled to the 
optical luminosity of a galaxy: slowly rising and low am- 
plitude for low-luminosity galaxies, high central gradient 
and high rotation velocities for high-luminos ity systems (e.g. 
iRubin et all Il985l ; iBurstein fc Rubinl H985h . However, the 
question of whether or not other optical properties influ- 
ence rot ation curves as we ll ha s resulted in inconsistent 
answers. IRubin et all (|l985h and IBurstein fc Rubinl (|l985l ) 
found no dependence on morphological type or on the shape 
of the light distribution (notably the bulge-to-disk ratio) 
and presented synthetic rotation curves depending solely 
on a galaxy's luminosit y. Th is idea was la t er ela borated by 
IPersic fc Saluccil f|l99lh and IPersic et all (|l996l ). who pre- 
sented a 'universal rotation curve', an analytic formula de- 
scribing the shape of a rotation curve which only depends 
on total luminosity. 

In contrast, several other studies suggested that ro- 
tation curve shape is not correlated with luminosity only, 
but tha t other parameters need to be taken into account 
as well. ICorradi fc Capacciolil (Il990l ) found that the shape 
of a rotation curve correlates with a galaxy's morpholog- 
ical type: early-type spirals with large bulges have rota- 
tion curves which rise more rapidly than galaxies of sim- 
ilar luminosity but with a less conce ntrated light distribu- 
tion . [Casertano^&^an Gorkoml (|l99lh showed that the outer 
shape of rotation curves is correlated with both the total 
luminosity and the shape of the light distribution, exem- 
plified by two luminous galaxies with highly concentrated 
light d is tribut ions which have declining rotation curves. 
iRoscoel (I1999T) showed that the universal rotation curve 
formalism of IPersic et al.l (| 19961) can be improved by in- 
cluding surface brightness as parameter influencing rota- 
tion curve shape. The dependence of rotation curve shape 
on the op tical characte ri stics was also c onfirmed in stud - 
ies by e.g. iBroeils! (tl992h. ISwatersI (Il999h. I Verheiienl (|200ll ). 
iMatthews fc GallagheT(|2002r ) and lSancisil (|2004f ). 

The systematics behind rotation curve shapes hold im- 
por tant clues on the st ructu re and evolution of (disk) galax- 
ies. IRubin et all (| 19851 ) and IBurstein fc Rubinl \l98$ ) inter- 
preted the lack of dependence on the light distribution as 
evidence that luminous matter plays a minor role dynami- 
cally, and that large amounts of dark matter must be present 
everywhere in disk galaxies. But if rotation curves are, in- 
stead, a multi-parameter family depending also on proper- 
ties such as morphological type, surface brightness, etc., then 
the conclusion must be that at least in some galaxies, the 
stars contribute significantly to the potential. The rotation 
curve vs. optical properties relations also provide a powerful 
benchmark for simulations of galaxy formation: any viable 
theory of galaxy formation must be able to reproduce real- 
istic rotation curves which match the other characteristics 
of the simulated galaxy. 

In order to obtain a better understanding of these is- 



sues, a systematic study of Hi rotation curves in spiral 
galaxies, covering a large range of luminosities, morpholog- 
ical types and surface brightnesses, is a crucial step. Al- 
though mu ch work has been done in this fie l d in recent 
years (e.g. Ide Blok et all Il996l: ISwaterj Il999l : ICote et all 
l200d : IVerheiienl l200ll : iGentile et al.l |2004 )~ most studies 
have focused on late-type and low-luminosity galaxies. 
Early-type disk galaxies , which generally contain less gas 
(|Roberts fc Hav ncs 1994; [Noordermeer et alj|2005l ). have re- 
ceived considerably less attention. One of the few studies 
so far aimed at a systematic investigation of Hi rotation 
cur ves ove r the fu ll range of morphological types was that 
bv iBroeilsl $mk. However, in his sample of 23 galaxies, 
only one was of morphological type earlier than Sb and only 
four had Knax > 250 km s" 1 . The only large-scale Hi survey 
dir ected spec i fically at SO and Sa galaxies was carried out 
bv Ivan Driell (|1987), but his study was severely hampered 
by the low signal-to-noise ratio of his data and his rotation 
curves were of rather poor quality compared to modern stan- 
dards. In the optical, little work has been done on early- type 
spiral galax ies either, si nce the early studies bv lRubin et al.l 
(1985) and I Kent! (1 19881 ). SO and Sa galaxies were th us also 
under-represented in the study bv IPersic et all (|l996h : their 
Universal Rotation Curve is based on over 1000 rotation 
curves of which only 2 are of type Sab or earlier. 

This paper is part of a larger study designed to fill 
this lack and to systematically investigate the relation be- 
tween dark and luminous matter in early-type disk galax- 
ies. These systems, lying at the high mass, high surface 
brightness end of the disk galaxy population, are ideal test 
cases to investigate what determines the shape of rotation 
curves. If the stars contribute significantly to the gravi- 
tational potentials of galaxies, it is in these galaxies that 
their influence will be most ea sily detected. In an earlier pa- 
per dNoordermeer et al.l l2005l. hereafter paper I), we have 
presented Hi observations for a sample of early-type (SO - 
Sab) disk g alaxies, and in an accompanying pa per to the 
present one l|Noordermeer fc van der Hulstll2006l . paper II) 
we present optical photometry and bulge-disk decomposi- 
tions. H ere, we use the data for a subset of 19 galaxies from 
iPaper il to derive their rotation curves and to study the de- 
pendence of their rotation curve shapes on the optical prop- 
erties. In two future publications, we will use the results to 
study the location of massive, early-type disk galaxies on the 
Tully-Fisher relation and to create detailed mass-models. 

The Hi data from I Paper il can be used to measure the 
rotation velocities of the gas out to large radii. In the cen- 
tral regions, however, the rotation curves can often not be 
measured from the 21cm o bservatio ns due to the presence 
of holes in the Hi disks (see lPaper J) . Furthermore, the spa- 
tial resolution of our Hi observations is usually insufficient 
to obtain detailed information on the shape of the rotation 
curves in the inner regions, where our velocity fields suffer 
from beam smearing. To overcome these difficulties, we use 
long-slit optical spectroscopy to measure the central rota- 
tion curves. In most galaxies, optical emission lines can be 
detected in the very inner regions, out to radii where reli- 
able rotation velocities can be determined from the Hi ve- 
locity fields. Moreover, due to the higher spatial resolution 
of the optical observations, the effects of beam smearing are 
strongly reduced. 

The remainder of this paper is structured as follows. 
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Table 1. Sample galaxies: basic data. (1) sample number; (2) UGC number; (3) alternative 
name; (4) morphological type; (5) distance; (6) and (7) absolute B- and R-band magnitudes 
(corrected for Galactic foreground extinction); (8) R-band central disk surface brightness (cor- 
rected for Galactic foreground extinction and incl ination effects) and (9) R- band disk scale 
length. Column (4) was taken from NED, (5) from lPaper il and (6) — (9) from Pap er ij . 



sample 


UGC 


alternative 


Type 


D 


M B 


Mr 




hn 


number 




name 




Mpc 


mag 


mag 


mag 
ai cscc^ 


kpc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1 


624 


NGC 338 


Sab 


65.1 


-20.83 


-22.25 


21.92 


5.8 


2 


2487 


NGC 1167 


SA0- 


67.4 


-21.88 


-23.24 


20.12 


8.0 


3 


2916 






Sab 


63.5 


-21.05 


-22.01 


20.99 


5.0 


1 


2953 


IC 356 


SA(s)ab pec 


15.1 


-21.22 


-22.54 


19.25 


4.1 


5 


3205 






Sab 


48.7 


-20.89 


-21.88 


19.59 


3.5 


6 


3546 


NGC 


2273 


SB(r)a 


27.3 


-20.02 


-21.35 


19.49 


2.8 


7 


3580 






SA(s)a pec: 


19.2 


-18.31 


-19.42 


21.58 


2.4 


8 


3993 






SO? 


61.9 


-20.19 


-21.35 


22.37 


5.5 


9 


4458 


NGC 


2599 


SAa 


64.2 


-21.38 


-22.61 


21.26 


8.6 


10 


4605 


NGC 


2654 


SBab: sp 


20.9 


-20.09+ 


_t 


_t 


_t 


11 


5253 


NGC 


2985 


(R')SA(rs)ab 


21.1 


-20.86 


-21.90 


21.32 


5.3 


12 


6786 


NGC 


3900 


SA(r)0+ 


25.9 


-19.94t 


-21.13 


19.30 


1.5 


13 


6787 


NGC 


3898 


SA(s)ab 


18.9 


-20.00 


-21.28 


20.49 


3.3 


11 


8699 


NGC 


5289 


(R)SABab: 


36.7 


-19.48 


-20.74 


22.24 


3.7 


15 


9133 


NGC 


5533 


SA(rs)ab 


54.3 


-21.22 


-22.62 


21.27 


9.1 


16 


11670 


NGC 


7013 


SA(r)0/a 


12.7 


-19.20 


-20.55 


19.58 


1.8 


17 


11852 






SBa? 


80.0 


-20.44 


-21.53 


20.74 


4.5 


18 


11914 


NGC 


7217 


(R)SA(r)ab 


14.9 


-20.27 


-21.35 


19.91 


2.7 


19 


12043 


NGC 


7286 


SO/a 


15.4 


-17.53 


-18.26 


19.90 


0.8 



t No data available in Paper II; Mb taken from LEDA. 



In section [2] the criteria which were used to select suitable 
galaxies from the parent sample of IPape^H are described. 
Section [3] describes the techniques that were used to derive 
the rotation curves from the Hi velocity fields and from the 
optical spectra. In section [4] the fitted orientation parame- 
ters and systemic velocities of our galaxies, as derived from 
different sources, are compared. In section [5] we briefly dis- 
cuss the occurrence of warps in the galaxies in our sample. In 
section H3 several aspects of the shape of our rotation curves 
are discussed, including an analysis of the correlations with 
optical properties and the applicability of the concept of a 
'Universal Rotation Curve' to our data. Finally, we briefly 
discuss our results and summarize the main conclusions in 
section [7] In the appendices, we present some additional 
material. A detailed description of the individual rotation 
curves is presented in appendix^] In appendix [B] we inter- 
pret the broad central velocity profiles which are present in 
some of our optical spectra. Appendix[C]grves the graphical 
representation of the rotation curves and various other data 
for the galaxies in our sample. 



be moving in regular circular orbits around the centre of the 
galaxy. Strongly interacting galaxies, or galaxies with oth- 
erwise distorted kinematics cannot be used. Strongly barred 
galaxies are excluded as well, because non-circular motions 
in the bar potential complicate the analysis of the data; 3) 
the inclination angle must be well constrained and prefer- 
ably lie between 40° an d 80°. 

Few galaxies from IPaper J satisfy all these conditions 
and a strict application of these criteria (especially the sec- 
ond one) would lead to a very small sample. We have there- 
fore relaxed the latter two selection criteria and included a 
number of galaxies with e.g. weak bars, mild kinematical dis- 
tortions or a more face- or edge-on orientation. The resulting 
sample consists of 19 galaxies; a few basic characteristics of 
the members are given in table [T] 

The galaxies in our sample have morphological types 
ranging from SO- to Sab and span two decades in optical 
luminosity ( — 17.5 > Mb > —22). The majority of galaxies 
in o ur samp le have high optical luminosity, with Mb < — 20. 
See lPaper il for a more detailed description of the properties 
of the galaxies in our sample. 



2 SAMPLE SELECTION 

The galaxies for the rotation curve study presented here 
were selec ted from the 68 galaxies with Hi observations pre- 
sented in IPaper j which were in turn selected from the 
WHISP survey (Westerbork survey of Hi in spiral and ir- 
regul ar galaxies; iKamphuis et al ] 1 19961 : 1 van der Hulst et al.l 
l200ll ). In order to be able to derive high quality Hi rotation 
curves, galaxies were selected on the basis of the following 
criteria: 1) the velocity field must be well resolved (> 5 - 10 
beams across) and defined over significant parts of the gas 
disks (i.e. not confined to small 'patches'); 2) the gas must 



3 OBSERVATIONS, DATA REDUCTION AND 
THE DERIVATION OF THE ROTATION 
CURVES 

As mentioned in the introduction, the rotation curves in 
this paper were derived from a combination of Hi synthesis 
observations and long-slit optical spectra. Below, we discuss 
the analysis of both components separately. 
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Table 2. Dynamical properties: (1) UGC number; (2) and (3) position of the dynamical centre; (4) heliocentric systemic velocity; 
(5) position angle (north through east) of major axis; (6) inclination angle; (7) maximum rotation velocity; (8) rotation velocity 
at 2.2 R-band disk scale lengths; (9) asymptotic rotation velocities at large radii; (10) total enclosed mass within last measured 
point and (11) rotation curve quality. 

UGC dynamical centre V sys PA i Vm ax U 2 . 2 h V" aS ymp M cnc quality 

RA (2000) Dec (2000) 





h 


m 


s 


O 




// 


km/s 







km/s 


km/s 


km/s 


M 




(1) 




(2) 






(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


624 


1 





36.0 


30 


40 


10 


4789 


288 


64T 


300 


300 


270 


5.2 ■ 10 11 


III 


2487 


3 


1 


42.7 


35 


12 


21 


4952 


250-256 


36 


390 


360 


330 


2.1 ■ 10 12 


I 


2916 


4 


2 


33.5 


71 


42 


19 


4537 


242 


42-50 


220 


210 


180 


2.8 • 10 11 


11 


2953 


4 


7 


46.8 


69 


48 


46 


892 


98-104 


50 


315 


315 


260 


1.1- 10 12 


I 


3205 


4 


56 


14.9 


30 


3 


8 


3586 


230-224 


67 


240 


230 


210 


4.3 ■ 10 11 


I 


3546 


6 


50 


8.6 


60 


50 


46 


1837 


56 


55 


260 


185 


190 


2.4 ■ 10 11 


II 


3580 


6 


55 


31.2 


69 


33 


54 


1203 


6-356 


63 


127 


100 


125 


9.0 • 10 10 


II 


3993 


7 


55 


44 


81 


55 


33 


4364 


220 


20 


300 


290 


250 


8.6 ■ 10 11 


11 


4458 


8 


32 


11.2 


22 


33 


36 


4756 


288-295 


25 


190 


280 


240 


7.8 ■ 10 11 


11 


4605 


8 


49 


11.1 


60 


13 


14 


1347 


240-247 


84-74 


225 


220* 


185 


2.8 • 10 11 


I 


5253 


9 


50 


22.2 


72 


16 


44 


1329 


356-340 


37 


255 


245 


210 


5.5 • 10 11 


II 


6786 


11 


19 


9.2 


27 


1 


15 


1795 


181-186 


68-64 


230 


_$ 


215 


3.1 ■ 10 11 


I 


6787 


11 


49 


15.3 


56 


5 


5 


1172 


107-118 


69-66 


270 


250 


250 


5.0 • 10 11 


I 


8699 


13 


45 


7.7 


41 


30 


19 


2516 


280 


73 


205 


190 


180 


1.9 • 10 11 


I 


9133 


14 


16 


7.7 


35 


20 


37 


3858 


24-45 


53 


300 


265 


225 


1.3 • 10 12 


I 


11670 


21 


3 


33.5 


29 


53 


50 


774 


336-330 


70-68 


190 


155 


160 


1.6 ■ 10 11 


II 


11852 


21 


55 


59.6 


27 


53 


55 


5843 


200-175 


50-60 


220 


210 


165 


5.9 • 10 11 


II 


11914 


22 


7 


52.3 


31 


21 


36 


951 


265-268 


31 


305 


300 


300 # 


1.9 • 10 11 


II 


12043 


22 


27 


50.4 


29 


5 


45 


1007 


97-92 


67 


93 


82 


90 


3.2 • 10 10 


I 



I Kinematical inclination poorly constrained by observations. Value copied from optical isophotal analysis. 

* No accurate photometry available due to edge-on orientation of optical disk; optical scale length is estimate only. 



" Galaxy does not have regular exponential disk; no optical scale length available. 

* Rotation curve extends out to 3.3 R disk scale lengths only and may converge to different velocity at larger radii. 



3.1 Hi rotation curves 

The Hi rotation curves were derived b y fitting tilted ring 
models (|Begem an 1987; Bcgcman 1989) to the observed ve- 
locity fields from lPaper A using the ROTCUR algorithm im- 
plemented in GIPSY (Gro ning en Image Processing System; 
IVogel aar fc Terlouwl 120011 ^1. In iPaperJ . we showed velocity 
fields at either full (« 15"), 30" or 60" resolution. Here, 
we fit tilted ring models to the velocity fields at all available 
resolutions. The higher-resolution velocity fields can be used 
for the inner regions, whereas the velocity fields at lower res- 
olution generally extend out to larger radii and can be used 
to obtain information about the rotation curves in the outer 
parts. Tilted rings were fitted to the entire velocity fields, 
but points near the major axis were given more weight than 
those near the minor axis by applying a |cos(a)| weighing 
scheme, with a the azimuthal angle, measured from the ma- 
jor axis in the plane of the galaxy. 

In all cases, the rotation curves were determined in four 
steps. In the first step, all parameters (i.e. systemic velocity 
Kys, dynamical centre position (x c ,y c ), position angle PA, 
inclination angle i and rotation velocity Kot) were left free 
for each ring. In general, the fitted systemic velocities and 
dynamical centre positions show little variation with radius, 
especially in the inner regions, and the average values were 
adopted as the global values for the galaxy. They are listed 
in tabled 



For two highly inclined galaxies, UGC 4605 and 8699, the stan- 
dard tilted ring method is not suitable and a modified analysis 
was applied (see individual notes in appendix [Aj . 



In the second step, the systemic velocity and dynami- 
cal centre were fixed for each ring at the values derived in 
the first step. The values for the position angle derived from 
this step are shown with the data points in the figures in 
appendix [C] The position angle is usually well-defined, but 
it often shows variations with radius as a result of warps 
in the gas disk. If a clear trend was visible, we fitted it by 
hand and used the fitted values for the next steps; otherwise 
we used the average of all rings. The adopted range of posi- 
tion angles, or the average value, for each galaxy is given in 
table[2] an d plotted as bold line in the figures in appendix [Cl 

In the third step, only the inclination and rotation ve- 
locity were left as free parameters for each ring. From this 
fit, the inclination angle was determined. This parameter is 
the most difficult one to constrain, because it is strongly 
coupled to the rot ation velocity, especially for in clinations 
lower than ~ 45 (|Begemanlll987l ; lBegemanlll989T ). The fit- 
ted values are shown in the figures in appendix [C] It is clear 
that for the more face-on galaxies, the uncertainties in the 
fitted inclinations are large. In practice, radial variations in 
inclination could only be detected for galaxies that are suf- 
ficiently inclined; for galaxies with i 45° only an average 
value could be determined. When necessary, the fitted incli- 
nation angles were also compare d to the values derived from 
the optical images (see lPaper ij) to make a more reliable es- 
timate. The range of inclination angles, or the average value, 
used for the next step is given in table[2] an< i plotted as bold 
line in the figures in appendix [Cl 

The uncertainty in the inclination Ai(r) was estimated 
by eye, based on the spread of the individual data points 
around the fitted values and the comparison between the 
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Table 3. Observational parameters for optical spectroscopic observations: (1) UGC number; (2) telescope 
used: Isaac Newton Telescope (INT) or William Herschel Telescope (WHT) on La Palma or NOAO 2.1m 
telescope on Kitt Peak (KP 2.1m); (3) observing dates; (4) total exposure time; (5) and (6) effective slit 
width and length; (7) position angle (north through east) of the slit on the sky and (8) - (12) line weights 
used in the stacking procedure. 



UGC telescope dates t cxp slit orientation line weights ui; 

width length PA [Nii] 6548 Ho [Nii] 6583 [Sn] 6 7i6 [Sn] 673 i 



(1) 


(2) 


(3) 


sec 
(4) 


// 

(5) 


(6) 


O 

(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


624 


INT 


26/1/01 


2400 


1.5 


3.3 


106 


0.5 


2.0 


1.0 


0.5 


0.5 


2487 


INT 


28/1/01 


7200 


1.0 


3.3 


70 


0.5 


1.0 


1.0 


0.5 


0.5 


2916 


INT 


27/1/01 


2400 


1.0 


3.0 


76 


0.0 


1.0 


1.0 


0.0 


0.0 


2953 


WHT 


2/1/00 


7200 


1.0 


4.0 


99 


0.5 


1.0 


1.0 


0.5 


0.5 


3205 


INT 


26/1/01 


6000 


1.5 


3.3 


17 


0.0 


1.0 


1.0 


0.0 


0.0 


3546 


INT 


26/1/01 


3600 


1.5 


3.3 


57 


0.0 


1.0 


2.0 


0.0 


0.0 


3580 


INT 


28/1/01 


3600 


1.0 


3.3 


5 


0.25 


1.0 


0.5 


0.5 


0.5 


3993 


INT 


26/1/01 


7200 


1.5 


3.0 


41 


0.25 


0.5 


1.0 


0.5 


0.5 


4458 


INT 


27/1/01 


3600 


1.0 


3.0 


100 


0.25 


1.0 


1.0 


0.25 


0.25 


4605 


INT 


28/1/01 


4800 


1.0 


3.3 


63 


0.0 


2.0 


1.0 


0.5 


0.5 


5253 


INT 


27/1/01 


2400 


1.0 


3.2 





0.1 


1.0 


1.0 


0.3 


0.3 


6786 


INT 


29/1/01 


2400 


1.0 


3.3 


2 


0.0 


l.ot 


l.Ot 


0.0 


0.0 


6787 


INT 


31/1/01 


6000 


1.0 


3.3 
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0.1 


1.0 


0.5 


0.25 


0.25 


8699 


INT 


22/5/01 


3600 


1.0 


3.6 


100 


0.25 


1.0 


1.0 


0.5 


0.5 


9133 


INT 


22/5/01 


7200 


1.0 


3.6 


27 


0.0 


1.0 


1.0 


0.0 


0.0 


11670 


INT 


22/5/01 


4200 


1.0 


3.6 


157 


0.25 


1.0 


1.0 


0.5 


0.5 


11852 


INT 


23/5/01 


4800 


1.0 


3.6 


15 


0.0 


1.0 


1.0 


0.0 


0.0 


11914 


INT 


23/5/01 


3150 


1.0 


3.6 


89 


0.25 


1.0 


1.0 


0.25 


0.25 


12043 


KP 2.1m 


8/12/01 


2400 


1.0 


5.2 


98 


0.25 


2.0 


0.5 


0.7 


0.0 



t Lines could not be stacked because of stellar absorption feature in Ha (see note in appendix I A|l . 



tilted ring inclination angles and optical ellipticity. In gen- 
eral, we let the uncertainty Ai increase with radius, in or- 
der to account for the possibility of undetected or misfitted 
warps in the outer gas disks of the galaxies. The adopted 
uncertainties in the inclination angle are shown with the 
shaded regions in the bottom middle panels in the figures in 
appendix [C] 

In the final step, we derived the rotation curves by doing 
a fit with all parameters fixed except the rotation velocity. 
To prevent the inclusion of erroneous points in the final ro- 
tation curves, outer tilted rings were not accepted if they 
only covered a small number of pixels in the velocity field, 
or if the outer parts of the velocity fields showed clear signs 
of non-circular or otherwise perturbed motions. The final 
Hi rotation curves are shown as square data points in the 
bottom right panels in the figures in appendix [Cl 



3.2 Rotation curves from optical spectra 

3.2.1 observations 

For most galaxies in the sample, long-slit spectra were taken 
with the IDS spectrograph on the Isaac Newton Telescope 
(INT) on La Palmfl For UGC 2953, a spectrum was ob- 
tained from the red arm of the ISIS spectrograph on the 
William Herschel Telescope, also on La PalmaQ The spec- 
trum of UGC 12043 was taken with the GoldCam spectro- 



graph, mounted on the NOAO 2.1m telescope on Kitt Peak, 
Arizonfl A summary of the observations is given in table [3] 

The slits of the spectrographs were aligned with the 
major axes of the galaxies. In some cases, the position angle 
of the slit on the sky was slightly different from the kine- 
matical position angle of the galaxy as derived from the Hi 
velocity field. In these cases, the rotation curves were later 
corrected for the effect of the misalignment. The bulges of 
the galaxies were usually bright enough to enable the slit 
to be positioned accurately over the centres using the TV 
camera in the focal plane. 

The spectral range of all observations was chosen such 
that each spectrum contains the redshifted lines of Ha (Ao = 
6562.80 A), [Nil] (6548.04 and 6583.46 A) and [Sn] (6716.44 
and 6730.81 A). The spectral resolution of the spectra taken 
on the INT is 1.0 and 1.4 A (FWHM) for slit widths of 1.0 
and 1.5" respectively, corresponding to a velocity resolution 
of approximately 45 and 65 km/s respectively. The spectrum 
for UGC 2953 has a spectral resolution of 0.9 A (~ 40 km/s), 
whereas the resolution of the spectrum for UGC 12043 is 
slightly worse at 2.0 A (~ 90 km/s). 

Total exposure times were broken up into single expo- 
sures of typically 20 minutes; the number of exposures for 
each galaxy was determined at the telescope, based on the 
strength of the emission lines in the first exposure. 



2 The Isaac Newton Telescope and William Herschel Telescope 
are operated on the island of La Palma by the Isaac Newton 
Group in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astrofisica de Canarias. 



3 The Kitt Peak 2.1m telescope is operated and IRAF is dis- 
tributed by the National Optical Astronomy Observatories, which 
are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Figure 1. Example spectra to illustrate the spectral stacking procedure. The top panel at each column shows a section of the original 
spectrum, centered around the Ha- (middle) and the 6583 and 6548 A [Nll]-lines (top and bottom respectively), after removal of the 
stellar continuum and sky lines. The middle panels show the stacked spectra, where the emission from all lines (including the [Sll]-lines) 
is added, with weights vii as given in table \3\ The bottom panels show the same spectra after binning along the spatial direction to 
~ 1" pixels. For each galaxy, the 3 spectra are shown on the same (logarithmic) intensity scale, to show the decrease in the noise levels 
between the subsequent steps. The arrows refer to specific features in the spectra and are discussed in the text. 



3.2.2 data reduction 

Standard data reduction steps were performed within the 
IRAF environment^. Readout bias was subtracted using the 
overscan region of the chips; any remaining structure was 
removed using special bias frames. The spectra were then 
flatfielded using Tungsten flatfields. Wavelength calibrations 
were performed using arc spectra from Copper, Neon and 
Argon lamps, taken before or after the galaxy spectra. The 
resulting wavelength solutions were used to map the spectra 
to a logarithmic wavelength grid, such that pixel shifts cor- 
respond to linear velocity shifts. The calibrations were also 
checked retrospectively by comparing the measured wave- 
len gths of a few strong nigh t-sky lines with the values given 
by lOsterbrock et alj (|l996l ); the systematic errors lie typi- 
cally in the range 0.05 - 0.10 A, corresponding to about 2.5 
- 5 km/s. 

Individual exposures were then combined and cosmic 
rays were rejected using a simple sigma-clipping criterion. 
The continuum emission of the galaxy and the night-sky 
emission lines were removed by fitting low-order polynomi- 
als along the spectral and spatial axes of the spectra respec- 
tively. 

In the final, cleaned spectra, the Ha line is usually the 
strongest line in the outer parts of the galaxies. In the central 
parts however, the 6583.46 A [Nil] line and the [Sll] lines are 
often stronger, presumably due to underlying stellar absorp- 
tion in Ha. Rather than first determining rotation curves for 
each line separately and then combining them into one sin- 
gle curve, we have chosen the reverse order. The parts of the 
spectrum around each of the 5 emission lines were shifted 
according to the difference in rest wavelength and stacked 
to create a single spectrum which contains emission from all 
lines. Each line was roughly weighted according to its rel- 
ative strength; the weights Wi are given in table [3] Before 
stacking, special care was taken to ensure that the different 
emission lines trace similar velocities, but we found no cases 
with significant differences. This procedure has the added 
advantage that the signal-to-noise ratio in the stacked spec- 
trum is higher than in the original one; emission that was 
too weak to be detected in each line individually could some- 
times be detected with sufficient significance in the stacked 
spectrum. 

In figure[T] we illustrate the stacking procedure for three 



representative galaxies. In the spectrum of UGC 624 (left), 
the Ha line is the strongest line along the entire slit, but 
adding the other lines leads to a slightly higher signal-to- 
noise ratio (see for example the location indicated with ar- 
rows 1), and thus improves the accuracy of the fitted ve- 
locities. For UGC 4605 (middle), the improvement is more 
significant. In the original spectrum, Ha is stronger in the 
outer parts, but the 6583.46 A [Nll]-line is stronger in the 
centre. In the stacked spectrum, velocities can be measured 
in both regions, as well as at locations where the signal in 
the individual lines was too weak to be fitted (arrows 2). For 
UGC 6786 (right), the stacking procedure does not work due 
to a strong stellar Ha absorption feature in the centre (arrow 
3). In this case, stacking the various lines causes the Ha ab- 
sorption feature to dilute the little emission that is present 
in the 6583.46 A [Nll]-line, and thus degrades, rather than 
improves, the quality of the data. In this case, we analysed 
both lines separately, and combined the resulting rotation 
curves afterwards (see also the note in appendix IA) . 

The final cleaned and stacked spectra are shown in the 
top middle panels in the figures in appendix [Cl 



3.2.3 derivation of the rotation curves 

From the stacked spectra, the radial velocity of the emitting 
gas was determined at each position along the slit by fitting 
Gaussian profiles along the wavelength direction. Before per- 
forming the fits, the spectra were binned in the spatial direc- 
tion to ~ 1 — 2" pixels to increase the signal-to-noise ratio of 
the data and to ensure that only one data point is fitted per 
resolution element. In some cases, parts of the spectra had 
such low-level emission that the signal-to-noise ratio was still 
too low in the binned spectra; for those regions, larger bin 
sizes were used. Spurious fits or fits with very large error- 
bars were discarded by hand. The fitted velocities are shown 
overplotted over the binned spectra in the top middle panels 
in the figures in appendix [C] They are also overplotted over 
a major-axis slice through the Hi data cube, shown in the 
top right panels in the same figures. 

The radial velocity curves for the approaching and re- 
ceding sides were then folded, using the centre of the op- 
tical continuum emission as central position. In two cases, 
UGC 3205 and 3580, the centres of symmetry of the emis- 
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sion lines appear shifted with respect to the location of 
the brightest continuum emission, in both cases by approx- 
imately one arcsecond (see the figures in appendix [C]). In 
these two cases, we determined by eye the central position 
which gave the largest degree of symmetry in the folded ro- 
tation curves. In the case of UGC 3580, the offset can easily 
be explained as a result of absorption of continuum emission 
by dust (see appendix lA")l ; for UGC 3205, the origin of the 
offset is unknown. 

The systemic velocity was determined by taking, at each 
radius, the midpoint of the velocities of the approaching and 
receding sides and taking the average of the resulting val- 
ues. This procedure maximizes the symmetry between the 
approaching and receding sides of the rotation curves over 
the full length of the spectra. In most cases, the systemic ve- 
locity thus derived is consistent with the value found from 
the Hi data (see also figure [2}. In some cases, small differ- 
ences were found; this happened mostly in galaxies which 
are kinematically lopsided, where an unambiguous determi- 
nation of the systemic velocity is difficult. In those cases, 
we closely inspected the optical spectrum and the Hi veloc- 
ity fields and determined interactively the systemic velocity 
which led to the smallest asymmetry in the final, combined 
optical and Hi rotation curves. 

Finally, the radial velocity relative to the systemic ve- 
locity was calculated for each point and, at radii where emis- 
sion was detected on both sides of the galaxy, the weighted 
average was determined. The final rotation curves were sub- 
sequently derived by correcting the average radial velocity 
curves for the inclination of the galaxy and for possible mis- 
alignments of the slit with the true major axis; the values for 
the inclination and position angle were taken from the re- 
sults of the tilted ring fits to the Hi velocity fields, described 
above. The resulting rotation velocities are shown with the 
filled circles in the bottom right panels in the figures in ap- 
pendix [C] 

3.2.4 optical beam smearing and other line-of-sight 
integration effects 

Close inspection of the optical spectra reveals that in many 
cases, the rotation curves rise so steeply in the centres of 
the galaxies that even in the optical spectra, the gradients 
are not fully resolved. Thus, although the optical spectra 
are a major improvement over the lower resolution of the Hi 
observations, they suffer from the optical equivalent of beam 
smearing as well and the fitted velocities in the central parts 
may still not represent the actual rotation velocities. 

Furthermore, many spectra have line profiles that are 
broadened even at positions several arcseconds away from 
the centres of the galaxies, where lack of resolution is not 
expected to play a major role anymore. These broadened 
profiles may be the result of line-of-sight integration effects 
through the disks and bulges of the galaxies. Again, the 
simple Gaussians which were fitted to these line profiles will 
not recover the true radial velocity at the projected radius 
and cannot be used for the final rotation curves. 

We have adopted a method similar t o the envelope- 
tracing (or terminal velocity) technique |Sancisi fc Alien! 
[Hzl ISofudl 19961 : iGarcfa-Ruiz e t~aTll2002l) to correct the in- 
ner points of the optical rotation curves which are affected 
by optical beam smearing and/or other line-of-sight inte- 



gration effects. We determined by eye the terminal veloc- 
ities of the affected line profiles, taking into account the 
instrumental velocity resolution. The effect of random mo- 
tions of the emitting gas clouds is ignored, as it is generally 
much smaller than the instrumental broadening of the pro- 
files (~10 vs. ~50 km/s). The radial velocities that were 
thus derived were then processed in the same way as the re- 
sults from the Gauss fits to derive the average inner rotation 
curve. The resulting rotation velocities are shown with the 
open circles in the figures in appendix [Cl 

Although the manually corrected rotation velocities are 
certainly a better approximation of the true velocities than 
the results of simple Gauss fits to the line profiles, there are 
many uncertainties, particularly regarding the detailed 3D 
distribution of the gas, that cannot be accounted for with 
the data used here. A more rigorous investigation of the 
kinematics in the central parts of the galaxies studied here 
would require even higher spatial resolution and preferably 
a fully 2D velocity field, i.e. either space-based or adaptive- 
optics assisted integral field spectroscopic observations. 

3.3 Final steps 

For the final rotation curves, the output from the tilted ring 
fits to the Hi velocity fields was compared to the derived op- 
tical rotation curves and it was determined which Hi data 
points were affected by beam smearing. Central Hi data 
points which lay significantly below the optical velocities 
were discarded. In almost all cases, the effect of beam smear- 
ing was limited to 1 - 2 Hi beam sizes from the centre, and 
only the inner two or three points of the Hi rotation curves 
had to be rejected. Only in highly inclined galaxies, such 
as UGC 4605 or 8699, do beam smearing and line-of-sight 
integration effects play a role at larger radii; these galaxies 
were treated individually to ensure that optimal corrections 
were applied (see appendix |X| . Outside the regions where 
the Hi observations are affected by beam smearing, the op- 
tical and Hi rotation curves generally agree to a high degree 
(<; 10 km/s). 

The remaining Hi data points were then combined with 
the optical data to produce the final rotation curves. Our 
final curves probe the rotation velocities over 2-3 decades 
of radii and enable us to measure small scale variations in 
the inner parts of the optical disks as well as the behaviour 
in the outer parts of the gas disks, many optical scale lengths 
away from the centre. 

The combined data points and their corresponding er- 
rors can be used, without further manipulation, to fit de- 
tailed mass models and to study the distribution of lumi- 
nous and dark matter in the galaxies; this will be done in 
a forthcoming publication. For the remainder of this paper, 
we are interested mainly in the global properties and shapes 
of the rotation curves. For this purpose, it is helpful to re- 
move the statistical fluctuations between the individual data 
points, especially those from the optical spectra. To do so, 
we fitted cubic splines through the data points, using the 
interactive fitting task CURFIT in IRAF. Individual data 
points from the rotation curves were weighted according to 
their errors (see below); points that were clearly offset from 
the main rotation curve were eliminated during the fits. The 
resulting curves are smooth but still follow the general be- 
haviour that underlies the individual data points; they will 
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be used in the remainder of the paper to study the shapes 
of the rotation curves and possible correlations with global 
properties of the galaxies (section [6]l . They are plotted as 
bold lines in the rotation curve panels in appendix [C] 

Finally, a few basic quantities are derived from the rota- 
tion curves. The rotation curves were classified on the basis 
of the quality and reliability of the data. Galaxies which are 
symmetric, show no signs of strong non-circular motion and 
have well-defined orientation angles are classed as category I. 
This class contains the following galaxies: UGC 2953, 3205, 
4605, 6786, 6787, 8699, 9133 and 12043. Category II contains 
galaxies with, for example, mild asymmetries, bar-induced 
streaming motions or signs of interactions or tidal distor- 
tions, as well as galaxies for which the orientation angles 
could not be constrained as accurately. The following galax- 
ies were classed as category II: UGC 2487 (Seyfert nucleus), 
2916 (interacting, lopsided), 3546 (strong bar and Seyfert 
nucleus), 3580 (lopsided), 3993 (inclination angle uncertain), 
4458 (possibly tidally disturbed), 5253 (tidally disturbed), 
11670 (large bar), 11852 (bar, kinematically disturbed) and 
11914 (inclination angle uncertain). UGC 624 was classi- 
fied as category III, because of the large-scale asymmetries 
present in the optical spectrum and particularly in the Hi ve- 
locity field. The rotation curve of this galaxy is of insufficient 
quality to be used for mass modelling. The classification of 
the rotation curves is listed in column (11) of table [2] 

We determined by eye the maximum and asymptotic 
rotation velocities, V ma x and Vasymp respectively. They are 
listed in table [2] and indicated with the horizontal arrows in 
the bottom right panels of the figures in appendix [Cl Simi- 
larly, the rotation velocity at 2.2 R-band disk scale lengths, 
V2.211 was determined. 

From the velocity at the outermost point of the rotation 
curve, the total enclosed mass is calculated as: 

M enc = = 2.325 • 10 5 ( fa* M , (1) 

G V km/s I kpc 

with V ou t the velocity at the last measured point and R ou t 
the corresponding radius. The resulting masses are listed in 
table[2]as well. In this calculation, it was implicitly assumed 
that the mass distribution interior to R ou t is spherical. If a 
significant fraction of the total mass is concentrated in a flat 
distribution, the value derived here is an upper limit. 



and line-of-sight integration effects. The contribution AV m 
is usually significant only in the optical rotation curves and 
in the inner parts of the Hi rotation curve, where only few 
points are available on the velocity field. At larger radii in 
the Hi rotation curves, where each tilted ring covers many 
data points in the velocity field, the measurement error is 
usually small (~1 — 2%). 

The second contribution AV nc comes from kinemati- 
cal asymmetries and non-circular motions in the galaxies. 
These were estimated by deriving rotation curves for the 
approaching and receding sides of the galaxies separately. 
Additional tilted ring models were fitted to the approaching 
and receding sides of the velocity fields and the resulting ro- 
tation curves were combined with the fitted velocities from 
the corresponding parts of the optical spectra. The result- 
ing rotation curves are shown with the crosses and plus-signs 
respectively in the bottom right panels in the figures in ap- 
pendix [C] The error in the rotation curve AVnc was then 
estimated as one fourth of the difference between t he rota- 
tion v elocities measured for each side separately (cf . ISwatersI 
1999). With this, rather ad hoc, assumption, the difference 
between the rotation velocity for each side separately and 
the average value represents a 2a deviation. Note that small- 
scale non-circular motions, or asymmetries perpendicular to 
the major-axis, are not accounted for in this estimate. 

The first two contributions to the rotation curve errors, 
AVm and AV nc were added quadratically, and are shown 
with the errorbars in the figures in appendix [Cl 

The third contribution to the rotation curve errors 
comes from the uncertainty in the orientation of the gas 
disks. The main contribution comes from the uncertainty Ai 
in the inclination angle, estimated as in section [3~T1 Errors in 
position angle are usually much smaller, and moreover, only 
contribute in second order to the rotation curve errors; in 
practice, they can be neglected compared to the uncertain- 
ties in inclination. The effect of the inclination errors on the 
rotation curves is derived as follows. The rotation velocities 
Kot in the rotation curve can be written as V TO t oc Kad/sin i, 
where V Ta d is the measured radial velocity from either the 
optical spectrum or the Hi velocity field. Thus, an error Ai 
in the inclination leads to an error AVi in the rotation ve- 
locity of 



3.4 Rotation curve errors 

Many factors can cause errors in the derived rotation ve- 
locities, both statistical and systematic. For a meaningful 
interpretation of the results, it is crucial to make a reliable 
estimate of all relevant uncertainties and much effort was 
therefore put into the identification and quantification of 
possible sources of errors. 

We account for three main contributions to the errors 
in the rotation curves. The first is simply the measurement 
error AVm. For the Hi data, this is given by the ROTCUR al- 
gorithm, based on the dispersion around the fitted tilted ring 
velocities; for the optical data, it is the fitted error on the 
profile centre, given by the Gaussian fitting routine. For the 
manually adapted velocities, the measurement errors were 
estimated by eye, based on the shape of the line profiles and 
the degree to which the data are degraded by beam smearing 



AVi = -^Ai md , (2) 



where Ai ra d is measured in radians. So, not only is it more 
difficult to derive the inclination accurately for near face-on 
galaxies, the resulting uncertainty in the rotation velocities 
due to a given error Ai becomes progressively larger as well. 

The derived errors AVi(r) are indicated with the shaded 
regions in the bottom right panels in the figures in ap- 
pendix [C] for clarity, they are drawn around the smoothed 
rotation curves, rather than around the individual data 
points. Note that these errors account not only for system- 
atic offsets of the rotation curves (due to a global misfit 
of the inclination), but also for the effect of undetected or 
misfitted warps, which would alter the shape of the rotation 
curves. 
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Figure 2. Comparison of global tilted ring (TR) parameters with values from other sources. Data points give offsets of central position 
(a), systemic velocity (b), inclination (c) and position angle (d) with respect to the values derived from the tilted ring analysis. For 
warped galaxies, the tilted ring values in the inner regions were used, such that they correspond to the sa me regions as probed by the 
optical observations. Filled squares represent the values from the optical isophotal analysis from lPaper III (a, c, d) or fr om the o ptical 
spectra (b). The crosses in panel b show the comparison with the systemic velocities derived from the global Hi profiles l|Paper fl ). The 
errorbars in panel c show the adopted uncertainties in the tilted ring inclinations. The standard deviations of the distributions are given 
in the bottom left corner of each panel. Cases where parameters from different sources were not independently derived (e.g. UGC 624) 
are not considered here. 



4 PARAMETER COMPARISON 

A proper derivation of a rotation curve depends crucially on 
the assumed orientation parameters and systemic velocity 
of the galaxy. It is instructive to compare the values which 
we assumed for the rotation curves here, derived from the 
tilted ring fits, with those obtained from other sources. 

Figure [2] shows that, in general, there is good agree- 
ment between the parameters from the tilted ring analysis 
and those derived from e.g. the optical isophotal analysis. In 
most cases, the dynamical and the isophotal centres coincide 
within a few arcseconds, well within one Hi beam. Larger off- 
sets are only observed in galaxies with strong dust features 
in the optical image (e.g. UGC 3580) or in highly inclined 
systems (e.g. UGC 8699), where extinction and line-of-sight 
integration effects complicate a proper determination of the 
central position. The observed offsets can be fully explained 
by observational effects, and no galaxies seem to have a true, 
physical offset between the dynamical and morphological 
centres. 

A comparison of the systemic velocities from different 
methods shows that they all agree within 5-10 km/s. In 
galaxies with well-resolved Hi velocity fields, however, the 
tilted ring systemic velocity is the preferred value, as it uses 



dynamical information from the entire gas disk. The other 
methods are expected to have larger intrinsic errors, so the 
dispersions given in the figure will probably come predomi- 
nantly from errors in those measurements; the typical error 
with which one can determine the systemic velocity using 
tilted ring models is probably of the order of 2 - 4 km/s. 

The disk orientation parameters derived from the op- 
tical isophotes usually agree within a few degrees with the 
tilted ring parameters. In particular, panel (c) shows that 
our assumed errors Ai on the inclinations are reasonable. 
Only 2 galaxies show an offset between the isophotal and 
kinematical inclination angles that is significantly larger 
than the assumed error. One is UGC 2916, which is interact- 
ing with a compani on galaxy . Close inspection of the fitted 
optical ellipticities IjPaper HT ) shows that the shape of the 
isophotes at intermediate radii is consistent with the kine- 
matical inclination; the outer isophotes are most likely dis- 
turbed by the tidal influence of the companion. The other 
case is UGC 6787, where the isophotal inclination is poorly 
constrained due to the influence of the dominant bulge. 
Again, since the tilted ring analysis uses dynamical infor- 
mation from the entire gas disks, it will generally give more 
accurate values for the inclination and position angles, so 
figure [2] mostly shows the uncertainties in the isophotal pa- 
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rameters. Note that the inclinations from LEDA have a large 
scatter around our values, with discrepancies up to 20°. 



5 WARPS 

It has been known for a long time that the outer parts of the 
gas disks of many spiral galaxies are n ot coplanar with the 
inner disk, b u t that they are 'warped' (jRoestad et al.lfl974l : 
ISancisilll976h . iBosmal i|l99ll ) reported th at at least 50% of 
all ga laxies are warped. More recently, iGarcfa-Ruiz et all 
(2002) studied 26 edge-on galaxies and found that all galax- 
ies with an Hi disk more extended than the stellar one are 
warped. 

Most galaxies in our sample are fairly face-on and warps 
are therefore seen less easily than in Garcia-Ruiz' galaxies. 
Nevertheless, we can infer the presence of warps from the 
tilted ring fits to the velocity fields. Inspection of the fig- 
ures in appendix [C] shows that the fitted inclination or posi- 
tion angles show significant radial variations in 14 of our 19 
galaxies. For three of the remaining five galaxies (UGC 624, 
3993 and 8699), the quality of the velocity fields is insuf- 
ficient to put strong constraints on the orientation of the 
gas disks and we cannot exclude the possibility that these 
systems are warped as well. Only two galaxies, UGC 3205 
and 3546, show little variation in the fitted orientation an- 
gles (the variations in the inner part of UGC 3205 can be 
attributed to bar-induced streaming motions) and seem to 
hav e no detectab le warp at all. 

claimed that warping tends to set in in the 
outer parts of the optical disk (around R25). Although many 
of the galaxies in our sample are consistent with having a 
flat gas disk within the optical radius, we also find a few 
counter-examples. The velocity fields of UGC 6786, 6787 and 
11852 show clear signs of warping in the inner parts. The 
first two systems are unbarred and the observed variations 
in their orientation parameters must be real. UGC 11852 
has a bar, but it is smaller than the Hi beam; the observed 
warping occurs at larger radii and must, again, be real. Note, 
however, that in all three cases, the inner warps are mild. 
Strong warps are only observed outside the bright optical 
disks, e.g. in UGC 9133 and 11852. 



6 ROTATION CURVE SHAPE 

The shape and amplitude of a galaxy's rotation curve are di- 
rectly related to the gravitational field in the midplane of its 
disk, and thus to the mass distribution of its main compo- 
nents. A systematic study of the shapes of rotation curves, 
and a comparison with the optical properties of the galaxies, 
can therefore yield important information on the distribu- 
tion and relative importance of dark matter in galaxies. In 
particular, many studies have addressed the correlation be- 
tween the distribution of luminous matter and the shape of 
a rotation curve (see the introduction for references). If the 
luminous matter plays a significant role in the dynamics of 
galaxies, the shape of a rotation curve must depend on the 
distribution of the luminous matter. On the other hand, if 
dark matter is dominant, such a correlation will be much 
weaker or completely absent. 

In figure O we show a compilation of all the rotation 



curves in our sample. In the left hand panels, the rota- 
tion curves are plotted on the same physical scale; in the 
right hand panels, all radii are sc aled with the R-band scale 
lengths of the stellar disks ffrom lPaper IJ) . Although there 
is a large variety in rotation curve shape among the galaxies 
in our sample, there are also some general features which 
can be recognised from this figure and from the individual 
rotation curves shown in appendix [Cl 

Almost all rotation curves in our sample rise extremely 
steeply in the central regions. In only one case (UGC 12043 
(#19)) do we see the 'standard' gradual solid-body-like rise 
of the rotation curve, before flattening out at about 3 disk 
scale lengths. In all other cases, the initial rise of the rota- 
tion curve is unresolved, even in the optical spectrum, and 
the rotation velocities rise from to J> 200 km/s within a 
few hundred parsecs (or, similarly, within a fraction of a 
disk scale length). In some cases (such as UGC 2953 (#4), 
3205 (#5) or 3580 (#7)), the steep central rise is followed 
by a more gentle increase before the maximum rotation ve- 
locity is reached; in other cases (e.g. UGC 4458 (#9), 5253 
(#11), 9133 (#15)), the rotation curve rises to its maximum 
immediately. 

At larger radii, many rotation curves show a marked 
decline. In several cases (for example UGC 2487 (#2), 
2916 (#3), 5253 (#11)), the rotation curves are more 
or less flat in the inner regions and the decline sets 
in quite suddenly around the edge of the optical disks 
(near R25)', this behaviour is s i milar to that in e.g. 
NGC 3992 (iBottema fc Verheiienl |2002| ) and NGC 5055 
l|Battaglia et al.l l2006h T But there are also cases where the 
rotation velocities start decreasing well inside the optical 
disk (e.g. UGC 2953 (#4), 9133 (#15)), or even right from 
the first point in the rotation curve (UGC 4458 (#9)). 

Although the total decline in the rotation curve can be 
large (more than 50% in the case of UGC 4458 (#9); ~ 25% 
for UGC 9133 (#15) and 11852 (#17)), all declining rota- 
tion curves appear to flatten out at large radii. No rotation 
curves are found with a fully Keplerian decline in the outer 
regions, indicating that we have not yet reached the point 
where the mass density becomes negligible. Thus, although 
the rotation curves of massive, early-type disk galaxies look 
remarkably different from those of later-type spiral galaxies 
at small and intermediate radii (with the latter generally 
lacking the ste ep rise in the centre and the decline at inter- 
mediate radii; Corradi & Capacciolil Il990l . ISpekkens et al.l 
120051 . ICatinella et~al]|2006h . they show the same 'flatness' 
in the outer regions, proving that they too must contain 
large quantities of dark matter. An interesting possible ex- 
ception is UGC 4458, whose rotation curve only flattens out 
in the very outer regions. Given the large uncertainties in 
the outer data points due to the face-on orientation, we can- 
not strictly rule out that this rotation curve keeps declining 
in Keplerian fashion. We will investigate this issue, and its 
implications for the dark matter content in this galaxy, in 
more detail in our subsequent paper on the mass modelling. 

It is worth mentioning that many galaxies show dis- 
tinct features in their rotation curves (e.g. UGC 6787 (#13), 
8699 (#14); see also the notes on individual cases in ap- 
pendix • Only few galaxies have smooth rotation curves 
without 'bumps' or 'wiggles' and the declines at intermedi- 
ate radii are rarely featureless and monotonous. Although 
these irregularities may sometimes be caused by e.g. noise 
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Figure 3. Compilation of rotation curves from all galaxies in the sample. The top panels show the rotation curves with V max > 250km/s, 
the bottom panels show those with Vmax < 250km/s. In the left hand panels, the rotation curves are plotted on the same physical scale; 
the insets show the inner 5 kpc of all curves. In the right hand panels, all radii are scaled with the R-band disk scale lengths. The curves 
are identified with the sample numbers from table [T] To limit confusion, a few curves are plotted using different linestyles: UGC 2916 
(#3, dot-dash), 3993 (#8, short dash), 6787 (#13, long dash) and 9133 (#15, dots). 



or non-circular motions of the gas, they can often be rec- 
ognized on both sides of the optical spectra or Hi velocity 
fields and must, in most cases, reflect small-scale features in 
the underlying mass distribution. In particular, we will show 
in our forthcoming publication on the mass models that the 
'wiggles' and the detailed shape of the drop-off in the rota- 
tion curves can, in some cases, be linked to features in the 
light or gas distributions and can be used to constrain the 
relative contributions of the luminous and dark matter in 
these galaxies. 

In some cases, such as UGC 2953 (#4), 3993 (#8) or 
11670 (#16), there are indications that the rotation curves 
start to rise again at the outer edges of the Hi disks. Whether 
this effect is real or an artefact in the data is hard to tell. The 
corresponding points in the Hi velocity fields were derived 
from low signal-to-noise ratio line profiles and have large un- 
certainties. Furthermore, we cannot exclude the possibility 
that the gas in the outer regions moves on non-circular or- 
bits, or that we have not determined the inclination of the 
orbits correctly. Follow-up observations at higher sensitivity 
are required to investigate this in more detail. 

6.1 Correlations with optical properties 

To investigate the dependence of rotation curve shape on the 
optical properties of the galaxies, we have ordered the rota- 
tion curves from our sample according to several parameters. 
In figures U and [S] we present a compilation of our rota- 
tion c urves in a similar fashion as lCasertano fc van Gorkoml 
(1991); the rotation curves are ordered according to the 



maximum rotation velocity Vm ax and the R-band disk scale 
length (figure [4]) or bulge-to-disk luminosity ratio (figure [5}. 
In figure [5] we have divided our galaxies into different sub- 
samples, according to several optical parameters, and plot 
the rotation curves for each subsample separately. 



6.1.1 inner rotation curves 

Early results by iRubin et al.l (|l985l ) showed that the inner 
shape of a rotation curve is coupled to a galaxy's lumi- 
nosity: bright galaxies have steeply rising rotation curves, 
whereas low-luminosity systems reach the maximum rota- 
tion velocity at relatively larger radii. This relation was 
later confirmed by several other studies (e.g. lBroeilslll992l ; 
iPersic et al.lll996l ; IVerheiienlll997l ; ISwaterslI 19991 ). The galax- 
ies with the lowest luminosity (and corresponding maximum 
velocity) in our sample (UGC 3580 (#7) and 12043 (#19)) 
follow this trend and have rotation curves which rise rela- 
tively slowly. In particular, UGC 12043 (#19) is the only 
galaxy in our sample which completely lacks the charac- 
teristic steep rise in the centre; instead, its rotation veloci- 
ties increase gradually, in solid-body fashion and only reach 
the maximum around 3 disk scale lengths. The remaining 
galaxies in our sample, however, seem to indicate that the 
systematic progression breaks down above a maximum rota- 
tion velocity of ~ 200 km/s (see figure [If. All galaxies with 
a rotation velocity larger than ~ 200 km/s have the char- 
acteristic steep rotation curve in the centre. Whether the 
rotation velocities continue to increase after this initial rise, 
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Figure 4. Compilation of rotation curves following Cascrtano fc van Gorkomllll99ll) . The origin of each rotation curve is placed according 
to the maximum rotation velocity V ma x and the R-band exponential scale length of the stellar disk. Dashed lines are used to indicate 
the origins for galaxies where the central rotation curve was not measured. The individual rotation curves are labelled using the sample 
numbers from table [T] and scaled in radius and velocity, as indicated with the small axes at the bottom right. The scale lengths for 
UGC 4605 and 6786 (#10 and #12) are estimates only, so their exact position in the figure is uncertain. 



or whether the maximum is reached in the very centre, does 
not seem to depend on the total luminosity of the galaxy. 

Instead, the shape of the rotation curve in the inner 
regions seems to depend more strongly on the concentra- 
tion of the stellar light distribution. This can be seen most 
clearly in figure [5] and in the bottom panels of figure [6] where 
the rotation curves are ordered according to the bulge-to- 
disk luminosity ratio and the more generic measure of light 
concentration Rso / R20 respectively. These figures show that 
the rotation curves of galaxies with faint bulges and a rel- 
atively diffuse stellar light distribution continue to rise af- 
ter the steep central part, and reach the maximum outside 



the bulge-dominated regions (e.g. UGC 2953 (#4) and 3205 
(#5)). On the other hand, the rotation curves of galaxies 
with highly concentrated light distributions rise to the maxi- 
mum immediately. This also explains why UGC 12043 (#19) 
has such a shallow central rotation curve: it has no bulge 
component at all. The only system with a small bulge which 
appears to reach its maximum rotation velocity at very small 
radii is UGC 3546 (#6), but this galaxy has a Seyfert nucleus 
which makes its central rotation velocities highly uncertain 
(see the errorbars in the figure in appendix[C]); it is well pos- 
sible that the rotation curve of this galaxy rises more slowly 
than we have derived here. 




BD ratio 

Figure 5. Compilation of rotation curves similar to figure [4] this time with the origin of each curve placed according to the maximum 
rotation velocity V max and the R-band bulge-to-disk luminosity ratio. UGC 4605 was not included in this figure, since no bulge-disk 
decomposition was available for this galaxy. 



Thus, our da ta a ppear at odds with the claim of 
iRubin et all (|l985h and iBurstein fc Rubinl l| 19851 ) that op- 
tical morphology does not influence the shape of a rotation 
curve and that large amounts of dark matter must be present 
at all radii. Our data indicate that at least the bulge stars 
have a strong influence on the central rotation curves, and 
suggest that they dominate the gra vitational potential in 
the inner regions (in agr eement with Corradi fc Capacciolil 



Il99d . IVerheiienlll997l and lSancisill2004T K 



6.1.2 outer rotation curves 

Inspection of figures [4] -[6] shows that the shape of the rota- 
tion curves in the outer parts is correlated with the luminos- 



ity of the galaxies: luminous galaxies are more likely to have 
a declining rot ation curve than low-luminosity syst ems (in 
agree ment with lCasertano fc van Gorkomlll99ll and lBroeilsl 
1992). This is shown in a more quantitative way in panel 
a) of figure [Jj where we plot the ratio of asymptotic and 
maximum rotation velocity as a function of total absolute 
magnitude. This figure shows that all early-type disk galax- 
ies with Mr < — 20 have at least a modest decline in their 
rotation curve. 

The strength of the decline shows, however, lit- 
tle dependence on the s hape of the light distribution. 
ICasertano fc van Gorkoml l|l99ll ) concluded that, in a sam- 
ple of galaxies of type Sb and later, the most strongly declin- 
ing rotation curves occur in systems with a compact light 
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Figure 6. Rotation curves ordered according to (from top to bottom) absolute magnitude of the entire galaxy, absolute magnitude of 
the stellar disk, absolute magnitude of the stellar bulge, central surface brightness of the stellar disk and compactness of the stellar 
light distribution (measured by the ratio of the effective radii i?go an d R20 wh ich contain respectively 80 and 20% of the light). All 
parameters are derived from the R-band images (see tables A3 and A4 in lPaper III) . The bins are chosen in order to distribute the rotation 
curves evenly over the panels. All rotation curves are labelled with the sample numbers from table [T] and scaled with the R-band disk 
scale lengths. Since no accurate photometry is available for UGC 4605 (#10), this galaxy is not included in these plots. Furthermore, 
UGC 12043 (#19) does not have a bulge component and is not included in any of the plots in the third row. 



distribution, where 'compact' in their terminology meant 
'small disk scale length'. Our data show that such a cor- 
relation does not exist for early-type disks, as our sample 
also contains a number of galaxies with large scale lengths 
which have falli ng rotation cu rves (see figure [3J). This is in 
agreement with iBroeils! l|l992t) , who also found a number of 
large galaxies with declining rotation curves. In fact, two of 
the galaxies in our sample with the most strongly declin- 
ing rotation curves (UGC 4458 (#9) and 9133 (#15)), have 



large scale lengths (8.6 and 9.1 kpc respectively). Panel b) in 
figure [7] shows that, if a trend with linear size of the galax- 
ies exists_j3 1 t_aU i jt_js_i^^ direction as observed 
bv lCasertano fc van Gorkoml (|l99lh : larger galaxies have on 
average more strongly declining rotation curves. 

No trend is seen when, instead of the disk scale length, 
we use the more generic parameter RS0/R20 to define the 
compactness of the stellar light distribution (bottom pan- 
els in figure [6] and panel c) in figure [7J| . Declining rotation 
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Figure 7. Strength of the decline in the rotation curves vs. total magnitude (a), scale length of the stellar disk (b) and ratio of the 
effective radii i?go and R20 (c) (all measured from the R-band images). The open symbol indicates UGC 11914; due to the small radius of 
the last measured point in the rotation curve, its asymptotic rotation velocity is poorly defined (see also appendix |Aj . Since no accurate 
photometry is available for UGC 4605, this galaxy is not included in these plots. 



curves are seen both in galaxies with a compact light distri- 
bution and in galaxies with a more diffuse stellar component 
(such as UGC 2487 (#2) or 2953 (#4)). Note, however, that 
according to this criterion, the galaxy with the strongest 
decline in its rotation curve is also the most concentrated: 
UGC 4458 (#9). 



6.2 The Universal Rotation Curve for early-type 
disk galaxies 

IPersic et al.l (|l996u claimed, based on a study of over 600 
optical rotation curves and a small number of Hi rotation 
curves, that the shape of a rotation curve is solely governed 
by the galaxy's luminosity and can be described by the fol- 
lowing simple formula: 



Kmc (x) = Uo P t {(0.72 + 0.44 log A) 



+1.6 e" 
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Here, x = R/R op t is the radius expressed in units of the op- 
tical radius i? op t, the radius encompassing 83% of the light. 
U pt is the rotation velocity at f? op t and A = Lb/L* b is the 
B-band luminosity of the galaxy scaled with L* . In princi- 
ple, Vopt can also be related to the luminosity via the Tully- 
Fisher relation, but since we are mostly interested in the 
shape of the rotation curve here, we empirically determine 
Vopt from our observed rotation curves. 

This Universal Rotation Curve (URC) has received con- 
siderable attention in the literature, as it implies (together 
with other scaling relations such as the Tully-Fisher rela- 
tion) a tight connection between the luminous and dark 
matter in galaxies and, as such, has i mportant consequences 
for the theory of galaxy formation jP alcanton et all 1 19971 : 
iHernandez fc Gilmorell 19981 : lElizondo et al.lll999h . However, 
from the observational point of view, no consensus has yet 
been reached concerning the general applicability of the con- 
cept of the URC to real galaxies. Although the URC seems 
to give a reasonable description of the general trends in ro- 
tation curve shapes, it was readily noted that individual 
rotation curves often show large deviations from the URC 
|Courteaulll997l : IVerheiienlll997l : IWillicklll999l : | Garrido et al.l 



2004) and that other parameters than luminosity must also 
influence a galaxy's rotation cur ve (e.g. surfa ce density, 
bulge-to-disk ratios, etc., see also iRosco e3 ll999l ). Our find- 
ings that the rotation curves of early-type disk galaxies have 
distinctly different shapes (steep central rise, decline at in- 
termediate radii) than those of later-type systems of similar 
luminosity, and that within our sample, the rotation curve 
shape is only weakly coupled to luminosity and rather de- 
pends on factors such as light concentration, bulge-disk ra- 
tio, etc., raise additional questions on the ability of the URC 
to describe rotation curves of all classes of disk galaxies. 

In figure [8l we compare the predicted rotation curves 
from equation|3]to our observed rotation curves. We did not 
measure the optical radius R op t for our galaxies, but use 
Rso, the radius c ontaining 80% of the light in the B-band 
(see table A3 in IPaper II ). as an approximation instead. 
All velocities are scaled with Vgo, the rotation velocity at 
Rso- UGC 4605 was omitted from the analysis, since no ac- 
curate photometric data were available. The most obvious 
result from figure [5] is that the URC completely fails to ac- 
count for the steep central rise in our rotation curves. In all 
galaxies, except the bulgeless system UGC 12043, the URC 
severely under-predicts the rotation velocities in the centre. 
For UGC 2953, it is hard to see the discrepancy in the figure 
due to the crowding of the points, but also in this case, the 
observed velocities inside 0.4 Rso lie far above the predicted 
curve. Thi s failure in the inne r regions is, however, not sur- 
prising, as lPersic et al.l (|l996l ) derived their relations specif- 
ically for disk-dominated galaxies and did not take bulges 
into account. 

In the outer parts, the agreement is good in some cases 
(e.g. UGC 2953, 3205, 9133), but there are also many galax- 
ies where the observed rotation curves have a markedly dif- 
ferent shape than the predictions from the URC. In particu- 
lar, the slope of the outer rotation curve seems to be poorly 
predicted by the URC, with the observed rotation curves 
flat instead of rising (e.g. UGC 6786, 11670, 12043), declin- 
ing instead of flat (UGC 11852) or not declining rapidly 
enough (UGC 2487). But some of these differences at larger 
radii may be related to the presence of the bulges as well. 
Not only do bulges influence the observed rotation velocities 
in the galaxies, they also change the shape of the predicted 
rotation curve by altering the total luminosity and the opti- 
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Figur e 8. Comparison of our rotation curves (black data points with error bars) with the Universal Rotation Curve from lPersic et all 
(1996) (bold black lines). Shaded regions give the uncertainties due to inclination errors. Red data points show the rotation curves after 
subtraction of the bulge component (shown with the dashed red lines; not for the bulgeless galaxy UGC 12043). The solid red line shows 
the URC for the disk component. All rotation curves are scaled with Rgo, the radius containing 80% of the total B-band light (from 
iPaper itf) and V80, the rotation velocity at Rgo- See text for details. 



cal radius; thus, they may cause discrepancies over the full 
length of the rotation curves. 

To investigate to what extent the difference between 
the observed rotation curves and the URC can be explained 
by the presence of the massive bulges in our galaxies, we 
have subtracted their contribution from the observed and 
predicted rotation curves; the results are shown in red in 
figure [8] The bulge contributions were taken from our mass 
models, which will be presented in a forthcoming publica- 
tion. The solid red lines show the predicted rotation curves 



from equation [3] now using the absol ute B-ban d magnitudes 
of the disk components (taken from lPaper III) . 

The correction for the bulge influence has indeed al- 
leviated some of the discrepancies, especially in the central 
parts; some galaxies show almost perfect agreement with the 
predicted rotation curves now (e.g. UGC 2953, 6786, 9133). 
However, even after the bulge corrections, many differences 
remain. Most rotation curves still rise more rapidly than 
predicted by the URC and reach a flat plateau where the 
URC is still rising (e.g. UGC 624, 3205, 3546, 3993, 4458, 
11670, 11852). Also in the outer regions, the observed slopes 
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of the rotation curves often still differ from the predicted 
ones, although the differences are generally smaller than in 
the original curves. 

In conclusion, the foregoing analysis confirms the de- 
pendence of rotation curve shape on morphological type, 
with bulge-dominated, early-type disk galaxies having dis- 
tinctly different rotation curves than late-type spirals of 
similar luminosity. Thus, a universal rotation curve that 
depends only on luminosity is inadequate to account for 
the observed diversity in rotation curv e shapes along the 
Hubble-sequence. Although the URC of iPersic et~ail (|l996l 1 
may have its virtue as an empirical tool to parameterise 
global trends of several properties of disk-dominated galax- 
ies (surface brightness, size, etc.) with luminosity and the 
reflection of those on the gravitational fields, additional pa- 
rameters are required to account for the presence of bulges 
in earlier-type disks. In addition to the deviations in the in- 
ner regions, we have shown that the detailed shape of the 
drop-off in the outer parts of our rotation curves is not well 
reproduced either. Clearly, real galaxies are more complex 
than the simple URC prescription suggests and other factors 
than luminosity must contribute to the detailed shape of a 
rotation curve as well. 



7 DISCUSSION AND CONCLUSIONS 

In this paper, we have derived rotation curves for a sam- 
ple of 19 early-type disk galaxies (SO - - Sab) spanning al- 
most 2 decades in optical luminosity. The majority of the 
galaxies are luminous, with Mb < —20. The rotation curves 
were derived from a combination of Hi synthesis observa- 
tions and long-slit optical spectroscopy of the ionised gas 
and probe the rotational velocities and mass distributions 
on scales ranging from 100 pc to 100 kpc. Almost all of the 
rotation curves share a number of properties, which appear 
to be typical for this type of galaxies. 

The rotation velocities generally rise rapidly and often 
reach values of 200 - 300 km/s (and up to 500 km/s for ex- 
treme cases such as UGC 4458) within a few hundred parsecs 
from the centres of the galaxies. After the initial steep rise, 
the rotation curves show a diversity in shapes. In some cases, 
the rotation velocities gradually increase further and reach 
the maximum at intermediate radii. In other cases, the ro- 
tation curves remain flat after the initial rise, or even start 
to decline immediately. This diversity in shape appears to 
be related to differences in the light distribution in these 
galaxies: galaxies with concentrated light distributions and 
luminous bulges generally reach the maximum rotation ve- 
locity at small radii, whereas galaxies with a more diffuse 
stellar component, or small bulges, generally have rotation 
curves which peak further out. 

At larger radii, most rotation curves decline, with the 
asymptotic rotation velocity typically 10 - 20% lower than 
the maximum. The strength of the decline is coupled to the 
luminosity of the galaxy, more luminous galaxies having on 
average m ore strongly declining rota ti on cu rves, in agree- 
ment with ICasertano fc van Gorkoinl (Il99l[). However, we 
cannot confirm another claim of Casertano Sz van Gorkoml . 
that declining rotation curves occur preferentially in galax- 
ies with a compact light distribution. By 'compact', these 
authors meant 'small disk scale length'. In agreement with 



iBroeilsl (|l992T ), our sample also contains a number of galax- 
ies with large scale lengths which have falling rotation 
curves. Interestingly, two recent studies l|Spekkens et al.l 
120051 ; [ Catinella et al.l 120061 ) showed that later-type galax- 
ies, even those with high optical luminosity, do in general 
not have declining rotation curves. In contrast, we find that 
declining rotation curves are a characteristic feature of mas- 
sive, ea rly-type disk galaxies. This seems to suggest that, al- 
though [CasertMo^^^rrGortora were correct to claim that 
the shape of the light distribution determines whether or 
not a galaxy has a declining rotation curve, the term 'com- 
pactness' must be interpreted as the presence of a light 
concentration (i.e. a bulge) in the centre, rather than a 
s mall scale in absolute term s. Note that the two galaxies 
in ICasertano fc van Gorkoinl with the most strongly declin- 
ing rotation c urves (NGC 2683 and NGC 3521) both have a 
sizable bulge l|Kentlll985T ). Within in our own sample, on the 
other hand, the strength of the decline in the rotation curves 
seems barely related to the bulge-to-disk luminosity ratio, 
or generic light concentration Rso/R20- Thus, although our 
data suggest that the bulge plays an important role in this 
issue, we conclude that it is not possible to extract a sin- 
gle parameter from the light distribution of a galaxy which 
uniquely determines whether or not it has a declining rota- 
tion curve. Rather, it must depend in a more subtle manner 
on the relative masses, and the details of the mass distri- 
butions, of the various luminous and non-luminous compo- 
nents in a galaxy. We will address this issue in more detail 
in our forthcoming publication on the mass models for our 
galaxies. 

It is important to note that we have not found any 
rotation curve which declines in Keplerian fashion. In fact, 
all rotation curves flatten out in the outer regions. Early- 
type disk galaxies, despite appearing dominated by luminous 
matter in the central parts, must also contain large amounts 
of dark matter to explain the shape of the rotation curves 
in the outer regions. 

Two low-luminosity galaxies, UGC 3580 and 12043, 
have a distinctly different kinematical structure than the 
other systems in our sample. The rotation curve of the for- 
mer does rise rapidly in the centre, but where in most other 
galaxies the rotation velocities decrease or remain constant 
at large radii, they continue to rise almost all the way till the 
outer point in the case of UGC 3580 (out to ~8.5 disk scale 
lengths). UGC 12043 completely lacks the steep central rise 
in the rotation curve; instead, its rotation velocities increase 
gradually, in solid-body fashion, before becoming constant 
outside approximately 3 disk scale lengths. The rotation 
curves of these two galaxies resemble those of typical late- 
type and dwarf gal axies which generally have slowly rising 
rotation curves too jBroeilslI 19921 ; ISwatersll 19991 ). UGC 3580 
and 12043 also have different optical mo rphologies than most 
galaxies in our sample (see iPaper III) and the conclusion 
seems justified that low-luminosity early-type disks form an 
entirely different class of galaxies. 

We have compared our rotation curves wit h the predic- 
tions from the Universal Rotation Curve from I Persic et aO 
(1996). Since their model contains only one parameter, the 
total luminosity, and does not include morphological type, 
their URC fails to account for the steep central rise in our 
observed rotation curves. These discrepancies are reduced, 
but not removed entirely, when we subtract the bulge influ- 
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ence from the rotation curves. Furthermore, there are also 
many differences between the observed and predicted rota- 
tion curves at larger radii. The concept of a Universal Rota- 
tion Curve which depends only on luminosity appears to be 
insufficient to account for the observed diversity in rotation 
curve shape; other factors must contribute to the detailed 
shape of a rotation curve as well. 

All in all, the results presented here show that rotation 
curves form a multi-parameter family. Although luminosity 
is clearly a major factor determining the shape of a rota- 
tion curve, other parameters are important too. In particu- 
lar, early-type disk galaxies have distinctly different rotation 
curves than their later-type counterparts, an effect which 
we have shown is mostly due to the presence of bulges in 
thes e systems. This is i n contrast with some previous claims 
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) that the shape of the rotation curves is determined by 



a galaxy's luminosity only and that the way the light is 
distributed has little influence. Our results contradict this 
and, in fact, indicate that above rotation velocities of about 
200 km/s, the total luminosity has little impact on the shape 
of the central rotation curve and that, instead, the shape of 
the stellar light distribution governs the dynamics in the 
inner parts. Our findings have important consequences for 
our understanding of the structure of galaxies. In partic- 
ular, our data strongly suggest that, at least in the central 
regions of the early-type galaxies presented in this study, the 
luminous matter dominates the gravitational potential, with 
dark matter only starting to play a role outside the bulge- 
dominated regions. We will investigate this issue in more 
detail in a forthcoming publication, where we construct de- 
tailed mass models for the galaxies in our sample and study 
the relation between luminous and dark matter in a more 
quantitative way. 
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APPENDIX A: NOTES ON INDIVIDUAL 
GALAXIES 

UGC 624 (NGC 338) has strongly lopsided kinematics, 
which makes it difficult to determine the systemic veloc- 
ity accurately. When the systemic velocity is left as a free 
parameter in the tilted ring fits, it shows a gradual decline 
of almost 50 km/s towards larger radii, in an attempt to 
symmetrize the rotation curve. The value of 4789 km/s min- 
imizes the asymmetries in the central parts, but even then 
there are differences between the approaching and receding 
sides of the optical spectrum. On the approaching side, the 
rotation curve rises rapidly to a more or less flat plateau 
at about 285 km/s, whereas on the receding side, the rota- 
tion curve rises more gradually to a peak of approximately 
310 km/s, after which it declines slowly. In the outer parts of 
the Hi disk (R <; 50 kpc), the differences are even more pro- 
nounced. On the receding side, the rotation curve declines 
gradually, whereas it starts to rise again on the approaching 
side. At the outermost point, the difference between the two 
halves amounts to almost 100 km/s. 

It seems plausible that the asymmetries in UGC 624 are 
caused by gravitational interaction with its neighbour, 
UGC 623. Note that the distrib ution of the neutral gas in 
UGC 624 is also asymmetric fcf. IPaper j) . 
The strong asymmetries in the velocity field lead to large 
residual velocities with respect to the model velocity field. 
Additionally, they make it difficult to accurately determine 
the inclination of the gas disk and we were forced to use 
the value from the optical isophotal analysis. This all re- 
sults in large uncertainties in the rotation curve, especially 
in the outer parts. This galaxy is therefore not suitable for 
a derivation of the dark matter properties and will not be 
used in our subsequent mass-modelling. 

UGC 2487 (NGC 1167) i s a giant SO galaxy (M B = 
-21.88, Di 5 = 54kpc; see IPaper III ) with an extended, 
highly regular gas disk. We can trace the Hi rotation curve 
out to radii of 80 kpc (10 R-band disk scale lengths) and al- 
though there is a small decline in the rotation velocities, they 
remain well above 300 km/s till the outermost point. The to- 
tal mass enclosed within the last measured point is M cnc = 
2.1-10 12 M Q , which makes UGC 2487 the most massive 
galaxy in our sample. The total enclosed mass is larger even 
than those in the giant Sc galaxies NGC 2916 and UGC 2885 
(|Rubin et all 1 19791 : iRoelfsema fc Alien! [19851 , note that in 
both papers a Hubble constant of Ho = 50 km s~ Mpc~ 
is assumed; their derived masses have to be divided by 1.5 
when using our value of 75 km s _1 Mpc -1 ); to our knowl- 
edge, it is the largest mass ever derived from a rotation 
curve. [Saglia fc Sanc isi (1988) list a number of other large 
disk galaxies with extremely high rotation velocities; some 
of those galaxies may be even more massive than UGC 2487, 
but since no spatially resolved rotation curves are available 
for these systems, no accurate values for the total masses can 
be derived. In any case, UGC 2487 seems member of a class 
of ex t remely massive disk g alaxies (see also iGiovanelli et all 
19861 : ICarignan et~allll997l ). with masses that rival those of 
the most massive e lliptical galaxies (e.g. iBertin et al1ll988l : 
iMinniti et ailll998l '). 

UGC 2487 is also classified as a Seyfert galaxy, ex- 
plaining the broad emiss ion lines in the nucleus (cf. 
iFilippenko fc Sargendll985l ). It has a central compact steep 



spectrum (CSS) radio source (e.g. ISanghera et all 1 19951 : 
iGiovannini et al]|200lj ). which is responsible for the Hi ab- 
sorption in the centre. Away from the bright nucleus, we 
detect some very faint emission in the optical spectrum. Al- 
though this emission seems to follow the general sense of 
rotation of the galaxy, the emission profiles are broad and 
do not have well-defined peaks. From these data alone, it 
is difficult to determine whether this faint emission traces 
regular rotation in the circumnuclear regions, or whether it 
is related to outflows from the active nucleus. Thus, this 
emission gives no useful information on the shape of the po- 
tential in the inner regions and we have decided not to use it 
in the derivation of the rotation curve. A small Hn region is 
detected 30" away from the centre on the approaching side; 
the emission from this region has regular line profiles and 
its velocity is consistent with the rotation velocities of the 
Hi at the corresponding location. 

UGC 2916 has a regular, symmetric rotation curve in its 
central regions. In the outer parts, however, a strong asym- 
metry is present between the approaching and the receding 
sides of the galaxy. At the approaching side, the rotation ve- 
locities show a strong increase outside the optical disk and 
reach a maximum of approximately 240 km/s at 20 kpc; fur- 
ther outwards the rotation velocities decline again. At the 
receding side, the rotation velocities do not rise at all outside 
the optical disk, but start to decline immediately. It seems 
likely that the lopsidedness in this galaxy is the result of the 
interaction with its nearby companion PGC 14370. 

UGC 2953 (IC 356) is by far the best resolved galaxy in 
our sample, with 117 independent data points in its rotation 
curve. The inner points, from the optical spectrum, sample 
the rotation velocities at intervals of ~80 pc or 0.02 R-band 
disk scale lengths, whereas the last measured point lies at 
a projected radius of 59 kpc (14 disk scale lengths). The 
central rise in the rotation curve is unresolved even in the 
optical spectrum; the rotation velocities rise to ~200 km/s 
within 1 arcsecond from the centre, and keep rising more 
gradually from there to a maximum of 310 km/s at i?~ 100" 
(~ 7.5 kpc, 2 disk scale lengths). In the outer regions, the 
rotation curve declines, with the asymptotic velocity about 
15% lower than the maximum. The strong variation in the 
fitted inclination angles around R = 600" is probably an 
artefact caused by streaming motions in the large spiral arm 
in the western parts of the galaxy; it was not judged to be 
real and was therefore not used in the final tilted ring fits. 

UGC 3205 has one of the most symmetric Hi disks of all 
galaxies i n our sa mple, both in its morphological appear- 
ance fsee lPaper J) and i n its kinem atics. The optical image 
is highly regular as well l|Paper III ). There is an almost per- 
fect symmetry in the velocity field, and the rotation curves 
for the approaching and receding sides of the galaxy are 
identical within the measurement errors. The residual veloc- 
ities are small too; the only significant residuals are detected 
in the regions where the bar causes non-circular streaming 
motions. It seems plausible that such streaming motions are 
also responsible for the apparent variation in inclination an- 
gle in the central parts; we find it unlikely that the fitted 
variation is real and have assumed a constant value for the 
final tilted ring fits. Outside the bar region, the residual ve- 
locities are of order 10 km/s or less everywhere, indicating 
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that the gas motion is highly regular and undisturbed. 
In the centre of the optical spectrum, a peculiar offset is 
detected between the centre of symmetry of the rotation 
curve and the centre of the continuum emission; the former 
is shifted by about one arcsecond (~ 250 pc) to the south- 
west with respect to the latter. It seems unlikely that this 
offset is caused by absorption of light by dust (as is most 
likely the case in UGC 3580), as the optical image is highly 
symmetric and shows no signs of dust extinction whatso- 
ever. No other peculiarities are seen in this galaxy at all and 
the origin of the offset remains unclear to us. For the deriva- 
tion of the optical rotation curve, we have used the centre of 
symmetry of the line emission, rather than the continuum 
centre, to fold the two halves onto each other. Given the 
flatness of the rotation curve and the high degree of symme- 
try at larger radii in the spectrum, a different choice for the 
dynamical centre would not have led to significantly differ- 
ent rotation velocities, except for the very inner points. 
The resulting rotation curve of this galaxy seems to lack, 
almost completely, the characteristic steep rise in the cen- 
tre which is observed in most other galaxies in our sample. 
Parts of the central regions are devoid of gas though, and 
we cannot trace the entire rise of the rotation curve. 

UGC 3546 (NGC 2273) is classified as a Seyfert galaxy 
l|Huchra et alj fl982). High-resolution observations with the 
VLA revealed a bright radio continuum source in the centre, 
consisting of two separate lobes separated by about .9" or 
«120 pc (jUlvestad fc Wilson! [l984 : iNaear et alJll999h . The 
optical spectrum shows strong nucle ar emission line s, the 
nature of which has been discussed bv lHo et all (|l997f ). The 
Ha emission in our optical spectru m is clearly extended (cf. 
|Poggelll989l : iMulchaev et~ai1ll996h . but the line profiles are 
irregular and broad out to a radius of about 10"; within this 
radius, it is difficult to disentangle the effects of quiescent 
rotation from possible in- and outflows from the nucleus. 
The uncertainties in the inner points of the rotation curve 
are therefore large, and the sharp peak and subsequent de- 
cline in the rotation curve may not be real. 
In contrast, the Hi velocity field displays smooth and reg- 
ular rotation throughout the entire gas disk of the galaxy. 
The central part of the velocity field appears somewhat dis- 
tort ed, but t his is an artefact caused by beam smearing (see 
also lPaper J) . This also explains the erratic behaviour of the 
fitted orientation angles in the inner regions; for the final 
tilted ring fits, we have used the position angle and incli- 
nation determined from larger radii, as indicated with the 
bold lines in the figures in appendix [C] No effects can be 
seen of non-circular motion in the bar, possibly due to its 
favourable orientation (perpendicular to the major axis). 

The centre of symmetry in the optical spectrum of 
UGC 3580 has an offset of approximately one arcsecond 
(~ 100 pc) with respect to the peak in the continuum emis- 
sion (indicated with the dashed line in the figure in ap- 
pendix [C]) . This difference is most likely caused by obscu- 
ration of the continuum emission by dust. As was noted in 
IPaper III , the optical image shows strong dust features in 
the central regions of this galaxy. Thus, a determination of 
the centre of the galaxy based on the peak in the light dis- 
tribution is problematic and it is not surprising that the 
dynamical centre is offset with respect to the isophotal one. 
Not only the optical appearance of this galaxy is peculiar, 



the kinematical structure is remarkably different from that 
of most galaxies in our sample as well. There is a marked 
asymmetry between the approaching and receding sides of 
the optical spectrum. Whether this reflects a true lopsided 
kinematics in the central regions or is a result of dust ex- 
tinction too, is difficult to determine from these data alone. 
The major axis slice through the Hi data cube is also asym- 
metric, but it appears that this asymmetry is peculiar to 
the major axis; it does not occur equally strong at different 
position angles and the rotation curves averaged over the 
full halves of the velocity field are only marginally different. 
Similarly to other galaxies in our sample, the central rise 
of the rotation curve of UGC 3580 is steep. But where in 
most of our galaxies the rotation curve becomes flat after 
the initial rise, or even starts to decline, it continues to rise 
gradually till twice the optical diameter (~8.5 R-band disk 
scale lengths) in this galaxy and only flattens out at the edge 
of the gas disk. This beha viour resembles that of typic al late- 
type and dwarf galaxies i|Broeilslll992l ; ISwatersi ri999) which 
generally have slowly rising rotation curves too. However, 
although UGC 3580 is one of the least luminous galaxies in 
our sample (Mb = —18.31) and correspondingly has one of 
the lowest maximum rotation velocities (Vmax = 127km/s), 
its luminosity is still too high to classify it as a dwarf 
galaxy. UGC 3580 thus seems to be a relatively luminous 
member of a class of low-luminosity early-type disk galax- 
ies which have distinctly different morphological and kine- 
matical features compared with their high- luminosity coun- 
terparts; other nice examples of this clas s are UGC 6742, 
UGC 12043 and UGC 12713 fsee lPaper J ). 

UGC 3993 is an SO galaxy with a regularly rotating gas 
disk. It resembles UGC 2487, but it is not as large nor as 
massive. Although the galaxy is quite face-on, the Hi veloc- 
ity field is of sufficient quality to determine the inclination 
with reasonable accuracy. Nevertheless, the uncertainties in 
the rotation velocities due to the inclination errors are large. 
In particular, we cannot exclude the possibility that the de- 
cline in the rotation curve of this galaxy is caused by a small 
warp in the outer parts towards a more face-on orientation. 

The rotation curve of UGC 4458 (NGC 2599) looks re- 
markably different from the canonical flat rotation curves 
normally observed in spiral galaxies. Instead, it rises very 
rapidly to an extreme peak velocity of 490 km/s, after 
which it shows a sharp decline of more than 50% before 
it asymptotically approaches a constant rotation velocity of 
~ 240 km/s. The central rise in the rotation curve is unre- 
solved even in our optical spectrum; it is possible that in 
the very inner regions, gas moves at even higher rotation 
velocities. Already, a rotation velocity of 490 km/s is unusu- 
ally high, and seems s urpassed only by that of UGC 12591 
|Giovanelli et alj|l986h . 

UGC 4458 is, however, close to face-on and the errors in the 
rotation velocities, caused by the uncertainties in the exact 
value of the inclination, are large. Thus, the peak velocity 
could be substantially lower if the galaxy were slightly more 
inclined in the inner parts, and the extreme decline in the 
rotation curve could partly be explained by a small warp to- 
wards a more face-on orientation in the outer parts. To fully 
explain the decline as a result of a warped gas disk, however, 
the inclination would have to decrease steadily from 31° in 
the centre to 15° in the outer parts (assuming a constant 
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rotation velocity of 400 km/s). Such a change in inclination 
angle is not detected in the tilted ring fits, and it seems un- 
likely that the decline in the rotation curve can be explained 
fully by warping of the gas disk. Meanwhile, we cannot ex- 
clude a small warp towards a higher inclination in the outer 
regions either. In particular, our data are also consistent 
with a continuing, Keplerian decline in the outer points of 
the rotation curve, and thus, with an absence of dark matter 
inside the radii probed by the Hi disk. We will investigate 
this issue in more detail in our subsequent paper on the mass 
modelling. 

The central parts of the gas disk of UGC 4605 (NGC 2654) 
are close to edge-on, with an estimated inclination angle 
of 84°. At larger radii, the gas disk warps towards a more 
face-on orientation. Due to the high inclination and the re- 
sulting line-of-sight integration effects, the central regions of 
the velocity field have a bias towards the systemic velocity 
and do not give an accurate r epresenta tion of the projected 
rotational velocities (see also IPaper J ), so the usual tilted 
ring analysis could not be used there. Within a radius of 
R — 30", this problem could be circumvented by using the 
optical spectrum; due to its higher spatial resolution, it suf- 
fers much less from projection effects. Between R = 30" and 
R — 110", where no optical information is available, we de- 
termined the rotation velocities by hand from the Hi data, 
using the same method that was normally used for the inner 
regions of the optical spectra (section |3.2[) . At each position 
along the major axis, the terminal velocity of the line pro- 
file was determined and assumed to represent the projected 
rotational velocity at the line of nodes; the average of the 
rotation velocities for the approaching and receding sides of 
the galaxy was then taken as the true rotation velocity at 
that radius. Outside R — 110", the gas disk becomes suf- 
ficiently less inclined that tilted ring fits could be applied 
to the velocity field; the rotation curve for R > 110" was 
therefore determined in the usual way. 

The edge-on orientation of this galaxy also has a positive 
effect, namely that the uncertainties in the inclination angle 
are small; the resulting errors in the rotation velocities are 
almost negligible. 

The final rotation curve is well defined and symmetric, ex- 
cept for the region around 120", where the rotation velocities 
on the approaching side are declining already while those on 
the receding side are still constant. Beyond this region, the 
rotation curve is symmetric and shows a strong decline on 
both sides. At larger radii (R <; 200"), the rotation curve 
flattens out at a level of approximately 185 km/s. 
Due to the edge-on orientation of this galaxy, no accurate 
optical photometry could be obtained and we were only able 
to obtain an estimate for the radial scale length of the stellar 
disk. But the rotation curve is relatively flat between 2 and 
3 times the estimated scale length, so the value for V2.2h is 
reasonably robust. Due to the lack of detailed information 
on the stellar mass distribution, however, this galaxy can 
not be used for mass modelling. 

The outer parts of the gas disk of UGC 5253 (NGC 2985) 
are dominated by a large spiral arm extending from the 
north of the galaxy. Although the gas in the arm is clearly 
rotating, it is impossible to determine the exact orientation 
of the arm and we have only fitted a tilted ring model to 
the inner parts. Even so, our rotation curve extends out to 



a radius of 49 kpc, or ~9 R-band disk scale lengths. 
The most interesting aspect of UGC 5253, however, is the 
strong m=0 component in the residual velocity field. Both in 
the inner (R <, 200") and in the outer (R <; 400") regions, 
the residual velocities are small, indicating that our fitted 
tilted ring model is an accurate description of the observed 
gas motions. Around a radius of 300" (~ 30kpc), however, 
a ring-like feature is detected in the residual field with an 
amplitude of —20 km/s. This feature has a high degree of 
symmetry with respect to the centre of the galaxy and the 
residual velocities are almost constant with position angle; 
only at the western end of the ring are the residuals slightly 
lower (~ — 10 km/s). In the rotation curve, the feature mani- 
fests itself as a marked asymmetry between the approaching 
and receding side. However, simple kinematical asymmetries 
can only explain residuals along the major axis of a galaxy; 
without additional radial motions, they must always vanish 
on the minor axis. The fact that the residual velocities in the 
ring are so symmetric and also present on the minor axis, 
argues against a simple explanation in terms of a kinemati- 
cal asymmetry. 

A more plausible explanation of the marked nr=0 component 
in the residual field is that it is induced by the gravitational 
perturbations fro m the large spiral arm in t he outer parts. 
It was shown by ISchoenmakers et al.l (l997) on theoretical 
grounds that an m=l perturbation in the gravitational po- 
tential leads to a strong m=0 term in the residual velocity 
field. The combination of the large m=l spiral arm in the gas 
distribution of UGC 5253 and the pronounced m=0 term in 
its residual field are thus a stro ng empirical confirmation o f 
their predictions. Furthermore, ISchoenmakers et al.l (|l997T l 
showed that it is, in principle, possible to use the amplitude 
of the m=0 term in the residual field to measure the strength 
of the perturbation on the potential. This analysis is, how- 
ever, beyond the scope of this paper and will be postponed 
to a later time. 

Alternatively, we could be seeing a vertical vibrational mode 
in the gas disk of UGC 5253, where the entire ring is mov- 
ing up wards (i.e. toward s us) with respect to the rest of the 
galaxy. ISellwoodl (ll996T l and lEdelsohn fc Elmegreenl (|l997f ) 
showed the results of N-body simulations which suggest that 
vertical vibrations can exist in the disks of galaxies, possibly 
triggered by a tidal interaction with a companion galaxy. It 
seems, however, questionable if such a mechanism could ex- 
plain the highly symmetric featu re we observe in UGC 5253 , 
and we judge the explanation of lSchoenmakers et al.l (|l997T ) 
more plausible. 

Note that the feature in the residual field coincides with a 
marked drop in the rotation curve. Both inside R = 200" 
and outside R = 400", the rotation curve is flat, but around 
R = 300", it suddenly drops from 245 to 210 km/s; the 
rate of the decline is consistent with pure Keplerian decay. 
Whether it is a coincidence that this drop occurs at similar 
radii as the feature in the residual field, or whether both 
effects are related, is unclear. 

A strong stellar Ha absorption feature is present in the 
central arcseconds of the optical spectrum of UGC 6786 
(NGC 3900); no Ha emission is detected in the inner parts. 
Emission is detected in the central parts in the 6583.46 A 
[Nil] line, but the standard procedure of stacking together 
the different emission lines in the optical spectrum cannot 



Rotation curves of early-type disk galaxies 23 



be used here. Instead, we analysed the Ha and [Nil] lines 
separately in this case, and combined the resulting rotation 
curves afterwards; at positions where emission was detected 
in both lines, the average velocity was calculated. The spec- 
trum shown in the figure in appendix [Cl was created by re- 
placing the inner part (15 pixels on either side of the centre) 
of the Ha spectrum by the corresponding region of the [Nil] 
line. Thus, it shows the extended Ha emission in the outer 
regions, together with the [Nil] emission in the nuclear re- 
gion. 

Note that the central emission in the [Nil] line is irregu- 
lar, with broad line profiles. This may be explained as a 
result of line-of-sight integration effects through the inne r 
regions of the massive bulge of UGC 6786 (see iPaper III) , 
but observations at higher spatial resolution and sensitivity 
are required to investigate this in more detail. At a radius 
of 5" on the approaching side, a strong emission feature is 
detected which has a velocity that lies more than 100 km/s 
closer to the systemic velocity than the emission at smaller 
radii. It seems unphysical that this emission traces regular 
rotation of gas in the plane of the galaxy and we did not 
include it in the combined rotation curve. 
The outer parts of the Hi disk are distorted as well, with 
two large spiral arms extending from the symmetric inner 
disk. It is impossible to determine the exact orientation of 
the gas in these arms and no tilted rings were fit outside 
a radius of 240". The residual velocity field shows a pecu- 
liar m=5 harmonic compo nent in the azimuth a l dire ction. 
According to the results of lSchoenmakers et al.l (|l997T ). this 
implies an m=6 perturbation in the gravitational potential 
of this galaxy. No obvious source of such a perturbation can 
be identified in the visible matter in this galaxy, so the per- 
turbation, if real, must be caused by the dark matter halo. 
Since UGC 6786 does not have a regular exponential stellar 
disk, no disk scale lengths are indicated in the figures in ap- 
pendix [C] 

The rotation curve of UGC 6787 (NGC 3898) is well re- 
solved and shows some characteristic 'wiggles' with an am- 
plitude of 30 - 50 km/s. The kinematics in the central parts 
are only barely resolved in the optical spectrum, and due to 
the high inclination angle and resulting line-of-sight integra- 
tion effects, the central line-profiles are strongly broadened. 
After the initial rise of the rotation curve to the peak velocity 
of 270 km/s, the rotation velocities drop to approximately 
220 km/s at a radius of 30" (~2.75 kpc), after which they 
gradually rise again to 250 km/s at R « 100" (~ 9 kpc). 
The rotation velocities then drop again to 220 km/s, after 
which they rise again to reach a more or less flat plateau at 
250 km/s. Although there are clear indications that the gas 
disk of UGC 6787 is warped, the locations of the 'wiggles' in 
the rotation curve do not coincide with the radii where the 
position angle and the inclination change and the variations 
in the rotation velocity seem real. This is further confirmed 
by the fact that the variations occur symmetrically at all 
position angles over the velocity field. 

The discrepancy between the kinematic incl ination an gle de- 
rived here and the optical inclination from IPaper ij can be 
explained by the domi nance of t he bulge in the optical im- 
age. As was noted in IPaper I H the optical image of this 
galaxy is dominated by the spheroidal bulge out to large 
radii, which makes it impossible to obtain a reliable estimate 



for the inclination from the isophotes. The kinematical in- 
clination derived here is free of such effects and thus reflects 
the true orientation of this galaxy more accurately. 

UGC 8699 (NGC 5289) is highly inclined and the central 
line profiles in the optical spectrum are severely broadened 
by the combined effect of beam smearing and line-of-sight in- 
tegration effects. Similar to UGC 6787, the rotation curve of 
this galaxy shows a distinct 'bump'. After an initial rise, the 
rotation curve reaches a peak velocity of 205 km/s at a ra- 
dius of approximately 8" (» 1.4 kpc). At R « 30" (» 5.5 kpc), 
the rotation velocities have fallen to 170 km/s, but unfortu- 
nately, no Ha emission is detected between 15 and 30" on 
either side of the galaxy, so the exact shape of the decline 
in the rotation curve cannot be recovered. Between R = 30" 
and R = 45", the rotation curve rises back to approximately 
200 km/s, symmetrically on both sides of the galaxy. At the 
approaching side of the optical spectrum, no emission is de- 
tected anymore beyond this radius; at the receding side, the 
spectrum indicates a small decline again in the rotation ve- 
locities, out to the last measured point at 54"(~9.7 kpc). 
In the Hi observations, UGC 8699 is poorly resolved along 
the minor axis. Comparison of the xv-slice through the Hi 
data cube with the optical spectrum shows that there is 
probably a central hole in the Hi disk, explaining the ab- 
sence of Hi emission at high rotation velocities close to the 
centre. Thus, this galaxy provides another good illustration 
of the use of the optical observations to resolve the shape 
of the inner rotation curve. Because of the poor resolution 
along the minor axis, standard tilted ring fits did not re- 
cover the rotation velocities accurately in this case. Instead, 
we used only points close to the major axis to determine the 
Hi rotation curve; points within 60° of the minor axis were 
discarded in the fits. In the outer parts, a small asymmetry is 
detected between the rotation velocities of the approaching 
and receding sides of the galaxy, with the former rotating 
about 10% slower than the latter. 

UGC 9133 (NGC 5533) has the most extended rotation 
curve in our sample, with the outermost point in the ro- 
tation curve lying at a projected radius of 103 kpc (> 11 
R-band disk scale lengths). With a rotation velocity at this 
radius of 225 km/s, the total enclosed mass is 1.3 ■ 10 12 M Q . 
Most gas at large radii lies in a giant spiral arm which ex- 
tends from the north-east side of the gas disk and is warped 
with respect to the inner parts of the galaxy. Although the 
arm is clearly rotating and the fitted position and inclina- 
tion angles are well-behaved, the lack of symmetry in the 
arm makes it difficult to exclude the possibility that the gas 
in the arm is not rotating at perfectly circular orbits. There- 
fore, care should be taken with the rotation velocities outside 
R « 200". Note, however, that we have assumed rather con- 
servative values for the uncertainty Ai in the inclinations, 
such that the corresponding uncertainties AVi in the rota- 
tion curve (indicated with the shaded area in the figures in 
appendix [Cj also include, at least partly, the uncertainties 
introduced by the orientation of the spiral arm. Note also 
that the residual velocities with respect to the tilted ring 
model are small, indicating that the non-circular motions 
are not dominant. 

The optical observations do not have sufficient resolution to 
trace the rise of the rotation velocities in the centre. Fur- 
thermore, a bright central component complicates the inter- 
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pretation of the spectrum in the inner few arcseconds (see 
also appendix [Bj. At larger radii, however, the spectrum is 
highly regular and symmetric. 

The rotation curve of UGC 9133 keeps declining almost to 
the outermost points in the Hi rotation curve. The asymp- 
totic velocity is about 25% lower than the maximum. Only 
at a radius of approximately 80 kpc does the rotation curve 
flatten out. 

UGC 11670 (NGC 7013) has a giant bar; the large 'wiggle' 
in the rotation curve around R = 50" may well be an arte- 
fact caused by streaming motions in this bar. The change in 
position angle which is detected at those radii is most likely 
also caused by non-circular motions and for the final rota- 
tion curve, we assumed that the position angle is constant 
at -24° in the inner regions. At larger radii, the rotation 
curve becomes smoother and shows a gradual decline out to 
a radius of about 300" (18 kpc), after which the rotation ve- 
locities seem to rise again. The upturn is only visible in the 
low-resolution data, but seems to occur on both sides of the 
galaxy, which strengthens the detection. We can, however, 
not exclude the possibility of in- and outflows in the outer 
parts of the gas disk of this galaxy; deeper observations are 
needed to verify if the gas in the outer parts is truly on cir- 
cular orbits and if the rise in the rotation curve is real. 
In the centre, a bright nuclear component complicates the 
interpretation of the optical spectrum (see appendix [B} and 
the inner points of the rotation curve may not be reliable. 

UGC 11852 has the most extended gas disk, relative to the 
optical size, of all galaxies in our sample. It allows us to trace 
the rotation curve to a radius of ~90 kpc, or ~21 R-band 
disk scale lengths. This galaxy is clearly warped: the fitted 
position angles from the tilted ring fits change by more than 
20° from the inner to the outer parts of the disk. The incli- 
nation changes as well and the galaxy becomes more edge-on 
in the outer parts. The apparent variation in the fitted incli- 
nations in the inner regions is probably an artefact caused 
by non-circular motions in and around the bar; it is not in- 
cluded in the final fit. The rotation curve can be seen to 
rise steeply in the inner regions, although the exact shape 
of the inner rotation curve cannot be recovered due to the 
strong central component present in the optical spectrum 
(see also appendix IB) . Outside the optical disk, the rotation 
curve shows a marked decline, but it seems to flatten out 
in the outer parts. The asymptotic rotation velocity is ap- 
proximately 25% lower than the maximum. Strong residual 
velocities (up to —50 km/s) are detected northwest of the 
centre of the galaxy. The velocities of the gas in this region 
are incompatible with regular rotation and indicate large- 
scale streaming motions. 

UGC 11914 (NGC 7217) has a highly regular gas disk; the 
rotation curves of the approaching and receding sides sepa- 
rately are almost indistinguishable. Unfortunately, however, 
the inclination is not tightly constrained and the resulting 
uncertainties in the rotation curve are large. 
Comparison of this galaxy with similar systems in our sam- 
ple (e.g. figures [3] -0 implies that it should have a declining 
rotation curve. The fact that we do not see such a decline 
is probably caused by the small radial extent of the Hi disk, 
which enabled us to measure the rotation curve out to 3.3 
R-band disk scale lengths only. In several other galaxies in 



our sample (e.g. UGC 2487, 5253), the decline in rotation 
velocities sets in around the edge of the optical disks, and 
it does not seem unreasonable to assume that the rotation 
curve of UGC 11914 declines at similar radii. 
Due to its relative proximity, UGC 11914 is one of the few 
galaxies where we fully resolve the inner rise of the rotation 
curve in the optical spectrum. We clearly see the gradual 
change in velocity from the approaching to the receding side 
over the central 30". Superposed over the 'normal' emission, 
we detec t the brigh t and unresolved central LINER compo- 
nent (cf. lKedll983h . The line profile of the LINER is very 
broad and implies velocities well in excess of the quiescent 
rotation velocities at larger radii (see also appendix iBjl. 

UGC 12043 (NGC 7286) is the least luminous galaxy in 
our rotation curve sample and the amplitude of the rotation 
curve is correspondingly low (Vm ax = 95 km/s). In contrast, 
the rotation curve is the second most extended, compared to 
the optical size, of all galaxies in our sample, with the out- 
ermost point lying at a projected radius almost 19 times 
larger than the R-band disk scale length. This galaxy is 
the only case in our sample with a slowly rising rotation 
curve. It completely lacks the steep central rise observed in 
all other galaxies in our sample and, instead, rises gradu- 
ally to the maximum which is only reached around 3 opti- 
cal scale lengths. This behaviour can be explained by the 
fact that UGC 12043 has no bulge and thus has a much 
smaller cent ral surfac e density than all other galaxies in our 
sample fsee iPaper III) . UGC 12043 is thus a member of a 
class of low-luminosity early-type disk galaxies which have 
distinctly different morphological and kinematical features 
than their high-luminosity counterparts; other nice exam- 
ples ofthisclass are UGC 3580, UGC 6742, and UGC 12713 
(see lPaper IT ). 

The interpretation of the optical spectrum is complicated 
by the low spectral resolution of the GoldCam spectrograph 
used for these observations. Although the rotation of the 
galaxy can clearly be detected, the line profiles are practi- 
cally unresolved and no corrections can be made for beam 
smearing or line-of-sight integration effects. A sudden jump 
occurs in the receding side of the optical spectrum, around 
R=20", the nature of which is unclear; no emission is de- 
tected at the corresponding radius on the approaching side. 
The Hi velocity field is highly symmetric, and the inferred 
rotation curves for the approaching and receding sides sep- 
arately are virtually indistinguishable, except for the outer 
parts of the low-resolution velocity field (not shown in the 
figure in appendix [Cj , where the approaching side exhibits 
a sudden rise in the rotation velocity. This rise is proba- 
bly caused by non-circular motions in the outer disk of the 
galaxy and most likely does not correspond to a real in- 
crease in the rotation velocity; it is not included in the final, 
azimuthally averaged rotation curve. 



APPENDIX B: BROAD CENTRAL VELOCITY 
PROFILES IN THE OPTICAL SPECTRA 

In about half of the galaxies in our sample, broad velocity 
profiles are present in the optical spectra of the central re- 
gions; the velocity amplitudes of those features often exceed 
the rotational velocities observed at larger radii. In some 
cases, especially in the Seyferts UGC 2487 and 3546, these 
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Table Bl. Central features in the optical spectra. (1) UGC number; (2) central velocity of 
line profile; (3) profile width, measured at 20% of the maximum flux; (4) derived rotational 
velocity; (5) estimated diameter; (6) total enclosed mass; (7) average density within D c ; 
(8) total R-band luminosity within D c and (9) R-band mass-to-light ratio within D c . 
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broad components may be explained as a result of nuclear 
activity and related gas flows. In the other cases, however, 
no strong nuclear activity is observed and the most natural 
explanation of the observed velocities is a nuclear disk or 
ring, possibly rotating around a massive nuclear star cluster 
or super-massive black hole. Similar rapidly rotating compo- 
nents have been observed in the center s of many other spi- 
ral galaxies (e.g. Carter fc Jenkinslll993|;lBertola et al.lll998l ; 
ISofue et al.lll99a Il999l; ISofue fe Rublj200ll and references 
therein; iMcDermid et al.ll2004l ; iFathil I2004T ). and may well 
be the spectro scopic counterparts to the nuclear structures 
imag ed by e.g. lCarollo et all (|l997l . ll99St ). 

In most cases, the central structures are spatially unre- 
solved in our spectra, but in a few cases (e.g. UGC 11670, 
11852) the resolution is sufficient to detect a velocity gra- 
dient. These gradients are always in the same direction as 
the sense of rotation of the gas at larger radii, which further 
strengthens our assumption that the central features in the 
spectra originate from regularly rotating gas. 

The data presented here lack the spatial resolution to 
unambiguously determine the nature of the central compo- 
nents in our spectra. A proper investigation would require 
sub-arcsecond resolution, both for the kinematic as for the 
photometric data, i.e. either space-based or adaptive-optics 
assisted observations. A nice example of t he potential of 
such o bservations was recently presented by I Atkinson et ah! 
(2005), who used imaging and long-slit spectroscopy from 
the Hubble Space Telescope to derive limits on the mass 
of the central black holes in NGC 1300 a nd 2748 (see also 
lHarms et ai]|l994l ; iFerrarese fc Fordlll999T l. 

With our data, we can only estimate the total mass 
in the central regions of our galaxies. The results of our 
crude analysis are summarized in table iBll The rotational 
velocities of the gas were derived from the width W c of the 
central line profiles, measured at 20% of the maximum flux, 
assuming that the central gas has the same orientation as 
the outer disk: V c = W c /(2smi). Note that the true ro- 
tational velocities may be larger or smaller if the central 
disks are tilted with respect to the main disk for which we 
derived the inclination angle, or if significant non-circular 
motions are present. The total extent of the gas in the cen- 
tral component was estimated by eye from the spectra di- 
rectly; table [Bll gives the estimated diameters D c . In most 
cases, the resolution of our data allows us to probe diame- 



ters of a few hundred parsecs or less; for the distant galaxy 
UGC 11852, the constraints are somewhat worse. Since in 
all cases, the structures are not or only marginally resolved, 
the diameters given in table IBll are highly uncertain and, 
strictly speaking, upper limits only. Finally, if we make the 
additional assumption that the central mass concentration 
has a spherically symmetric shape, the total enclosed mass 
can be estimated using equation [T] 

Table [Bll clearly shows that the inferred central masses 
and densities are high, especially in the nearby galaxies 
where we have better constraints on the total extent of the 
rapidly rotating gas. Central densities of order 10 3 M pc -3 
and higher can almost certainly not be explained by the 
normal stellar components. In column (8) of table IBll we 
give the total R-band luminosities wit hin D c , a s measured 
from the optical images presented in iPaper III It is clear 
that the observed luminosities are too small to account for 
the derived dynamical masses; the local mass-to-light ratios 
are much larger than what is expected for the surrounding 
bulge material. 

Our inferred masses are also a few orders of magni- 
tudes larger than thos e of the most massive stellar clus- 
ters known to date (e.g. Mengel et al. 2002; Mar aston et al.l 



(e.g.lMe 

120041 ; IWalcher et~aT] 120051 ). so it is unlikely that the high 
rotation velocities are caused by unresolved central con- 
centrations of stars. For UGC 2953 , 5253, 11670 and 
11914, HST images a re available (e.g. ICarollo et all |2002| ; 
iHunt fc Malkanll2004h and we could explicitly verify that no 
bright and compact sources of light are hidden in our own, 
lower-resolution images. 

In conclusion, the central components in our spectra 
seem to be a strong indication for the presence of super- 
massive black holes in at least a fraction of our galaxies. 
Sub-arcsecond observations, however, are required to obtain 
more detailed knowledge on the spatial extent and orienta- 
tion of the rapidly rotating gas and to provide conclusive 
evidence for the presence and mass of the black holes. 



APPENDIX C: ROTATION CURVES 

On the following pages, we present the rotation curves and 
several other related quantities for all galaxies in the sample. 
For each galaxy, we show a figure consisting of the following 
panels: 
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Left hand column: velocity fields 

Top: grayscale and co ntour rep resentation of the observed 
Hi velocity field from IPaper j Darker shading and white 
contours indicate the receding side of the galaxy. Contours 
are spaced at intervals of 25 km/s; the thick contour indi- 
cates the systemic velocity from tabled 
Middle: model velocity field, based on the tilted ring fits to 
the velocity field shown above. Contours and grayscales are 
the same as in the observed velocity field. 
Bottom: Residual velocity field, produced by subtracting the 
model velocity field from the observed one. The grayscales 
range from -40 to +40 km/s (white to black). Contours are 
spaced at intervals of 25 km/s; the thick contour indicates 
zero residual velocity. 

The cross in each panel indicates the dynamical centre, as 
given in table [2] 

Top row: x-v diagrams 

Middle: cleaned and stacked optical spectrum (see sec- 
tion [332] f° r a description of how these were produced). 
Right: major axis position-velocity slice through the Hi data 
cube. 

The position angles on the sky of both plots are indicated 
in the top left corner of the panels. Contours in both pan- 
els are at -1.5 and -3 (dotted) and 1.5, 3, 6, 12, . . . times the 
rms noise in the respective datasets. The dashed horizontal 
and vertical lines denote the systemic velocity and the cen- 
tre of the galaxy respectively. The fitted rotation velocities 
are overplotted. For each of the two sides of the plots, ap- 
proaching and receding, we plot the velocities derived from 
fits to that side only. Plotting symbols are as follows: 

- blue/black squares show velocities derived from the final 
tilted ring fits to the Hi velocity fields; 

- red/dark gray bullets show velocities fitted to the optical 
line profiles. 

- orange/light gray bullets show the central locations in 
the optical spectra which were affected by 'optical beam- 
smearing' and where the rotational velocities were de- 
termined by eye, rather than by Gaussian fits (see sec- 
tion [323J. 



the top show the corresponding radii in kiloparsecs. 
Plotting symbols and linestyles are as follows: 

- filled blue/black squares show velocities derived from the 
final tilted ring fits to the Hi velocity fields; the position and 
inclination angles assumed for these final fits are indicated 
with the bold lines in the middle panels of the bottom row. 

- open blue/black squares show velocities that were derived 
from Hi data at lower resolution than those shown in the left 
column and top row. 

- red/dark gray bullets show the velocities from the optical 
spectra. 

- orange/light gray bullets show optical velocities that were 
manually adapted to correct for beam-smearing and line- 
of-sight integration effects in the optical spectra (see sec- 
tion 

- crosses and plus-signs indicate the rotation curves for the 
approaching and receding sides respectively. 

- errorbars are a combination of fitting errors and differ- 
ences between the approaching and receding sides (see sec- 
tion [Mil. 

- bold lines show the smoothed rotation curves, derived 
from cubic spline fits through the individual data points 
(see section [3731 for details). 

- dashed lines indicate regions where no optical emission 
was detected. 

- shaded regions show the uncertainties in the rotational 
velocities due to the adopted uncertainties in the inclination; 
for clarity, they are drawn around the smoothed rotation 
curve, rather than around the individual data points. 

- vertical arrows at the bottom show the radii correspond- 
ing to one (left) a nd two ( right) disk scale lengths of the 
R-band image fsee lPaper lit) ; horizontal arrows at the right 
show the derived maximum and asymptotic rotation veloci- 
ties and the rotation velocity at 2.2 R-band disk scale lengths 
(see table [2]|. 



Bottom row, middle panels: Orientation parameters 
from the tilted ring fits to the Hi velocity fields. 
Top: data points show the fitted position angles (north 
through east) from fits with position and inclination angle 
left free. The bold lines give the values that were used for 
the final fits to derive the rota tion curves. The arrows give 
the values derived in IPaper III from the outer isophotes of 
the optical image. 

Bottom: data points show the fitted inclination angles from 
fits with inclination angle left free and position angle fixed 
at the values shown with the bold line in the top panels. 
The bold lines give the values for the inclination angle that 
were used for the final fits to derive the rotation curves; the 
shaded regions show the adopted uncertainties. The arrows 
give the values derived from the outer isophotes of the opti- 
cal image. 

Bottom row, right hand panels: Rotation curves. 
Left: inner regions. 
Right: full rotation curves. 

The axes at the bottom show radii in arcseconds, those at 
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